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Abstract 

A new method of including the complete NLO QCD corrections to hard processes 
in the LO parton-shower Monte Carlo (PSMC) is presented. This method, called 
KrkNLO, requires the use of parton distribution functions in a dedicated Monte 
Carlo (MC) factorization scheme, which is also discussed in this paper. In the 
future, it may simplify introduction of the NNLO corrections to hard processes and 
the NLO corrections to PSMC. Details of the method and numerical examples of 
its practical implementation as well as comparisons with other calculations, such 
as MCFM, MC@NLO, POWHEG, for single Z/7*-boson production at the LHC are 
presented. 
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1 Introduction 


Higher-order perturbative corrections in Quantum Chromodynamics (QCD), important 
for the LHC data analysis, are calculated order by order in the strong coupling, a s , 
while some numerically important ones, related to soft and collinear singularities, can be 
resumed to the infinite order. The most valuable, albeit technically difficult, way of QCD 
resummation is in form of a Monte Carlo (MC) event generator [I]. It is widely recognized 
that the most promising way of getting high precision QCD calculation for hadron collider 
data analysis is a common framework of fixed-order calculations combined with a parton 
shower Monte Carlo (PSMC). The pioneering work, in which the complete first-order QCD 
corrections to the hard process of heavy boson production in hadron-hadron collision were 
combined with PSMC, was that of Ref. |2|. Shortly later another interesting variant was 
proposed in Ref. [31. Presently both methods are available for many processes, see Ref. |4|. 

It is worth to mention that partial efforts in this direction were done earlier in the 
course of development of the most popular PSMCs, like for instance in Ref. (5] or |6| and 
in many other works. In these earlier attempts, the distributions generated in PSMC 
were improved using a tree-level exact matrix element (ME) of QCD, while virtual cor¬ 
rections were neglected or taken in the leading-logarithmic approximation |7|. It was also 
known for a long time, that an ad hoc approach, in which PSMC differential distribu¬ 
tions were corrected using the exact ME and the overall normalization was corrected by 
hand to fixed-order next-to-leading (NLO) integrated cross section, was providing distri¬ 
butions in a quite good agreement with experimental data. It may be therefore a little 
bit surprising that it took two decades to work out a systematic method of combining 
the NLO-corrections to the hard process, known from the early 1980’s, see for instance 
Ref. |[8], with the leading-order (LO) PSMC, also dating from the early 1980’s. Apart 
from the lack of interest in more precise QCD calculations due to poor data quality, 
main reasons for this much delayed development can be seen from problems addressed 
in Refs. [2,3]. Namely, any such a method requires a very good NLO-level analytic un¬ 
derstanding/control of distributions from PSMC and, either NLO-level complete phase 
space coverage for the hard process or a practical methodology of correcting for the lack 
of it. Luckily, a new wave of developments of the LO parton showers, see Refs. |9 12 


has lead to modernized PSMCs, better suited for merging/matching with the fixed-order 
QCD calculations, in particular with better or even complete coverage of the hard process 
phase space. 

It is now obvious that the next challenge on the way to even higher-precision pertur¬ 
bative QCD calculations needed until the end of the LHC era two decades from now, is to 
combine the fully exclusive NNLO corrections to the hard process and the NLO parton 
shower. The fixed-order NNLO corrections to many processes are already well established, 
see for instance Refs. 13,14, but the NLO PSMC needed for such a progress is still not 
available, except of feasibility studies summarized in Ref. fl5|. Interesting partial solu¬ 


tion of combining the NNLO-corrected hard process with the LO parton shower can be 
found irQrefs. 16-18 . The present work is relevant for the above future developments it 


1 In these methods only certain selected important distributions are upgraded to the NNLO level. 


2 











the sense that it presents a simplified method of correcting the hard process to the NLO 
level in combination with the LO parton shower (PS). In other words, it offers a simpler 
alternative to the MCONLO and POWHEG methods of Refs. |2,3|, which may hopefully 
pave the way to the NNLO hard process combined with NLO PSMC. 

The new method described here, nicknamed as KrkNLO, was already proposed in 
Ref. 19 , where its first numerical implementation was done on top of the dedicated 


toy model PSMC and was limited to the gluonstrahlung subset of the NLO corrections. 
Later on it was tested numerically in a more detail in Refs. 15 2 
the KrkNLO method is implemented within the standard PSMC Sherpa 2.0.0 


In the present work 
A 


11 


pilot study of KrkNLO implementation outside PS MC, using MC event encoded in the 
event records produced by Herwig++ 2.7.0 (9, 21, 221 and Sherpa 2.0.0 was also done. 
Let us stress, however, that the overall simplifications of the KrkNLO method comes not 
completely for free, as it requires to use parton distribution functions (PDFs) in a special 
Monte Carlo (MC) factorization scheme (obtained, however, easily from reprocessing the 
MS PDFs), and it is required that the basic LO PSMC provides for the NLO-complete 
coverage of the hard process phase space (this is also not a problem for all modern PSMCs). 
Our method is simpler to implement in the case of PSMC with an ordering based on the 
transverse momentum kr or a q 2 variable of Ref. 

Seymour work 


24 


23 , inspired by the classic Catani- 


However, it can be also easily implemented on top of PSMC that 
uses the angular ordering - without the need of the so-called truncated showers required 
in the POWHEG method, see Refs. 20,25 for more discussion on that. 


The main advantage of the KrkNLO method is a simplification of the NLO corrections 
due to the use of PDFs in the MC factorization scheme. The implementation of the entire 
NLO corrections with the help of a single multiplicative simple weight on top of the LO 
distribution is a quite unique feature of the KrkNLO method. 

Numerical studies presented here will be extended in the future publications to a 
wider range of distributions, energies, implementation variants, including comparisons 
with experimental data. 

The paper is organized as follows. In Section 2 we introduce the kinematics and the 
phase space parametrization for the considered process. In Section 3 we describe in detail 
the KrkNLO method. Sections 4 and 5 contain some numerical results of the KrkNLO 
implementation: Section 4 from the fixed-order NLO cross-checks while Section 5 from 
the NLO+PSMC comparisons with the MCFM, MCONLO and POWHEG programs for the 
main Z-boson observables. Section 6 summarizes the paper. In Appendix A we add some 
details on the first gluon emission in the backward evolution in PSMC. 


2 Kinematics and phase space parametrization 

In the present work we are going to concentrate on the Drell-Yan process, specifically 
production and decay of the heavy bosorj^] Z/q* in proton-proton collisions. At the 

2 For brevity, in what follows, we shall often speak about the Z boson only, but in all cases we really 
mean Z/q* —> l + l~. 
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Figure 1: Kinematics of Z-boson production in the qq channel. 

leading order (LO), qq —* Z is the only partonic subprocess that contributes. At the 
next-to-leading (NLO) level, the real correction qq —y Zg and the virtual correction to 
qq —>■ Z contributes - to be referred to collectively as qq channel. The remaining NLO 
contributions, qg —>• Zq and qg —>■ Zq, are tree-level only - to be called the qg channel. 


2.1 Single emission 

Fig. 0 illustrates part of the notation that will be used throughout the paper. The 
diagram shows the real correction to the Z-boson production in the qq channel with the 
gluon four-momentum denoted by k±. The four-momenta of the incoming forward and 
backward partons, p \f and p±s are related to the four momenta of the incoming protons, 
Pp and Pb, as follows 


POF — XfPf , 


POB — XbPb , 


( 2 . 1 ) 


where 

= 1,0 t ,1), Pb = ^(1,0 t ,-1) . (2.2) 

The invariant masses of the incoming parton pair prior to gluon emission and that of 
the produced Z boson are denoted by sq and s — s\, respectively. Their raticj^] 




(2.3) 


can be related to light-cone variables of the emitted gluon, as seen from the kinematics 
of the emitted gluon expressed in terms of the light-cone Sudakov variables 


2fci ■ Po B 
— - — 

s 0 

zi = 1 ~ «i - /3i , 


2 kf _2k 1 - p 0p _ 2 k x 

? Pi 

Sq Sq Sq 


kf T = aifaso , 


1 au 
Vl = 2 111 ft’ 


q?i + /3i < 1, ai, /31 > 0, 

s o — ^Po f Pq b i 


(2.4) 


3 The subscript in k\ and si is kept to underline that, in the context of PSMC, there is more parton 
emissions further away from the hard process qqZ vertex. 
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Figure 2: The Sudakov logarithmic plane with lines marking qf = const, (thick green) and 
Z\ = const, (thin blue). The dashed (red) line marks IR cut-off on k\r = const.. The shaded 
(yellow color) triangle marks the part of the phase space where ~ mp dominates. 


spanned on the four-momenta p\ F and p\ B of the incoming partons, prior to the gluon 
emission. As seen in Eqs. (2.4), the variables ay, Pi are simply the fractions of the 


light-cone components of the gluon four-momentum to the centre-of-mass (CM) energy 
of the incoming partons (prior to the gluon emission). The ratio z\, the gluon transverse 
momentum k\x and the gluon rapidity y\ are related to the ay, Pi variables as well. 

In PSMC, the (eikonal) phase space measure of the emitted gluon 


d 3 ki 

2fcf 


k 2 


7r dcj )i da\ d/3i 

2 2n ay P 1 


(2.5) 


19 , while 


is always split one way or another into two parts belonging to the quark and antiquark 
emitters. A sharp division along the y\ = 0 angular boundary was used in Ref. 
in modern PSMCs, see Refs. (9. TT| 12], a more gentle division, introduced in Refs. 

is exploited: 


and inspired by the Catani-Seymour work 24 


26,27 


da i d/3i 
«i A 


(m F +m B ) = 


daidfli 


da\df3i 


Pi(oii + Pi) au(aq + /3i) 


m F = 




Oil Pi 


m B = 


Pi 


a i + Pi 


( 2 , 6 ) 


where m FjB are the so-called soft partition functions. What is important for our method¬ 
ology in the following is that the sum of the two parts attributed to two emitters/showers 
reproduces the original 1-particle phase space without any gaps and (or well-controlled) 


5 





































overlaps. Clearly the rrip +mg = 1 property of the functions used in Eq. (2.6) takes care 
of that. 

Once the above (overlapping) separation of the emission phase space into two part is 
applied, a different evolution variable of PSMC in each of them is defined 


— s o(«i +/3i)/3i, ql B — so(aq + 0i)aq, (2-7) 

instead of a common one, like fcf. = aq0iSo or rapidity. 


is easily re-parametrized 


Using the above evolution variable and z±, the single-emission (eikonal) phase space (2.5 


doixdfii dq\ p dz x t dq\ B dz x 

oti/3i qf F 1 — zi 


+ 


dl 1 


Zl 


( 2 . 8 ) 


The relation between the old and new variables are illustrated graphically in Fig. [2] The 
transformation back to the Sudakov variables is different for each part: 


0i = 


aq = 


(Qi f /sq) 
1 - Z\ 
(sL/so) 


Oil = 1 — Zi *■ 
01 = 1 - Z! ~ 


1 ~ Z! 

The upper phase-space limit aq + 0i < 1 transforms into 


(<fi F /so) 

1 - U 

(gi B / g °) 

1 - A 


(2.9) 


Zi > 0 


and 


£ s o, 


( 2 . 10 ) 


while the positivity conditions, aq > 0 and 0i > 0, enforce the IR-boundary cut-offs 

(1 - Sl ) 2 > or (1-zr) 2 >^h, (2.11) 

SO s 0 

for the two parts, correspondingly. Also, in most of the phase space region populated 
according to the ~ m F factor we may approximate qf ~ k'f T . Similarly q\ ~ kf T in the 
~ m B phase-space sector. The above kinematical limits are also shown on the logarithmic 
Sudakov plane in Fig. [2] for the ~ mp sector. NB. The IR cut-off kf T > k'^ min marked in 
this figure translates into a slightly stronger cut-off on 1 — z±, easily calculable. The same 
kinematical limits for one emission are also illustrated directly in terms of the 1 — z± and 
q\ p variables in Fig. [3j including also the second emission for the purpose of the following 
discussion. 

The essential ingredient of the 1-particle phase space reorganization towards PSMC 
is introduction of the on-shell “effective beams” pi F and pi B , such that p\ F + p± B = 
Po F + Po B ~ k\. Their definition is not unique. For example, for the ~ mg branch one 
may choose 


Pi B = (! - e )Po B , Pi F =Po F - ki + ep 0B , e = 0i/(l-ai). (2.12) 
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hard 


Figure 3: The illustration of the kinematic boundaries in the forward evolution (FEV) algorithm 
with n — 2 particles. The parton no. 1 is emitted as a first one, within the biggest shaded 
triangular (blue) area defined by (1 — z\) 2 > qf/so and qf > Ql The parton no. 2 is generated 
as a smaller triangular area marked with dashed (red) line according to (1 — z 2 ) 2 > q 2 /s\ 
and q 2 > Qi- Third parton is not generated, but its would-be-allowed space is marked by 
the smallest triangle below s 2 = s . The vertical line marking the upper phase-space boundary 
q 2 < Q 2 in D(Q 2 ,x ) of Eq. (3.24) is also marked. 


The same for the ~ branch modulo the obvious variables interchange. All that in the 
rest frame of P 0 = po F + po B ■ However, in practice of a typical PSMC with the backward 
evolution, all four-momenta are reconstructed starting from qf 1 and z\ variables in the 
rest frame of the effective beams P\ = p± F +pi B , which are constructed in the first place. 

Technical details of the construction of effective beams are not so important for our 
analysis. Just as an illustration let us define explicitly one possible construction, which 
was introduced in Ref. [23]. In the rest frame of the hard process and the effective beams 
Pi — Pi F + Pi B = (Vi, 0,0,0) one may construct all four-momenta - starting from the 
(qi F , Zi, (pi) set, then translating it into («i,/3i) and using (for m F )\ 

K = O^lPlp + 0*Pl B + <?T, 1, Qt, 1 = <X*PlS, 

Po B = ~~~Pi B > Pop = Pi + k i- Po B , (2.13) 

* = 1-Pl-Z! = 1 Pi 

Pi + Zl 1 Pi + Z\ Zl 

Finally, knowing P 0 = Pi — k\ one may transform all newly constructed four-momenta 
k±, po F , Po B to the rest frame of Pq. 

Altogether, the complete reorganization of the 1-real emission phase space from a 
simple form based on the Sudakov variables to an equivalent parametrization using the 
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variables of PSMC, based on the backward evolution algorithm (applying the Catani- 
Seymour soft-partition factor mp), keeping track of the kinematical limits, and defining 
x — Xi — XqZi, s — Si — sx i, looks as follows: 


dcr\ F — dxiD(fi 2 F ,xi) 


dfiidai 


“i+/3i<l 


/3i(«i + Pi) 


P(zi) dcr 0 (sxoZi) 


1 so !—V^fTso 

= J dxi J J dzi dcTo(sxi) D(h 2 f ,xi) 


(2.14) 


= fdx J ^ J dzi 9 { 1 - Zl ) 2 s£/Zl > q 2 P q q{zi) da 0 (sx ) d(ji 2 f , 


0 (J 2 . 

’mi 


XQ 


where q\ = ql p , P(z) = (1 — z)P qq {z) and P qq (z) is the DGLAP 28 splitting function, 
ft is important to see that the full phase coverage requires integration over q 2 p to extend 
above the effective mass squared s = S 2 of the LO hard process. Since q p ~ k\ Tl it means 
that the transverse momentum above s is included in the above phase space. 

In Section [5] presenting numerical results, it will be commented more on what kind of 
evolution variable is chosen in the parton shower generation of Sherpa and Herwig++. 


2.2 Multiple emissions 

Although the above 1-emission kinematics is enough for most of our prescription for the 
NLO-correcting of the hard process, for better understanding of the role of the underlying 
multi-emission PSMC, it is useful to extend the 1-emission treatment of the kinematics 
to two and more emissions in the initial-state parton showers. 


BEV labeling 



Figure 4: Labelling in the reconstruction of the four-momenta. 


The important technical point is the choice of the numbering (labelling) of the emitted 
particles. In Fig. [3] we use the numbering in the emission chain (ladder) starting at the 
incoming hadron and ending next to the hard process, that is the labelling of the forward 
evolution (FEV) algorithm. From now on we switch to the backward evolution (BEV) 
labelling which starts next to the hard process and ends at the incoming hadron, see Fig. [4] 
for illustration. 










Let us summarize briefly on the effective beam technique , for simplicity limiting its 
description to a single tree F of emissions (shower). The sequence of the effective emitter 
beams (pi F ,pi B ) is defined starting from the hard process, with the four-momentum P t , 
such that 

i 

Pi = Pip + Pig, Pi+1 = Pi + k, Pf = S/ ( U ZiJ , 

3 = 1 

and they are used to span the four-momentum of the emitted gluon 

, ~ ~ s _ , _ kipi,, ~ kjpi„ 

k = otiPi F + Pip iB + kn, Oii = _ _ , Pi = . . , 

Pi F Pig PipPig 

introducing the Sudakov variables cq, Pi. These Sudakov variables are related to the 
evolution variable qf p and the light-cone variable Zi of PSMC as follows: 


P* 


f[Zi 

s(l - 


oti 


1 — Zi — Pi, Zi = YI Z F S = sx = sx 0 . 

3 = 1 


(2.15) 


Finally, in the recursive backward reconstruction of the four-momenta starting from 
the hard process, one employs the Sudakov-like decomposition of the emitted parton in 
terms of the emitters after the emissiorj^j 


K = a* 


= Y 


+ Qtai 


p{i-l)p+^i-l) B 

a* n* _ 1 A 

? Pi -1 

— oti 1 — Oif Zi 


(2.16) 


All exact kinematical limits (including ordering in the evolution variable) are represented 
by the following inequalities: 


(1 ~ Zi) 2 >q 2 iZi/s and q- +l < qf < ,§(1 - zrf/Z,. (2.17) 


In particular, the kinematical limits for the first emission in the backward-evolution (BEV) 
labelling are 

<?max,l = s(l - ZlY/Zl, Z?** - 1 “ 9o/V^- 

The above kinematical limits in terms of the BEV variables qf and Zi look more compli¬ 
cated than in the FEV scenario, although they represent exactly the same phase space 
region, and are illustrated graphically in Fig. [5] 

As already underlined, in terms of the BEV variables, the phase space gets apparently 
widened after each emission, for instance q 2 > s is already available for the first emission 
and, due to lowering of the IR boundary by the \jz\ factor, more phase space is available 
for the second emission. This phenomenon, important for the full coverage of the phase 


‘Generalizing Eq. (2.13). 
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Figure 5: The illustration of the kinematic boundaries in the BEV algorithm with n = 2 
particles. The parton no. 1 is emitted within the semi-triangular shaded (blue) area defined 
by (1 — Zi) 2 > qf/(s 0 /zi) and s > q\ > q 2 . The parton no. 2 is generated within the second 
semi-trapezoid area marked by the dashed (red) line according to (1 — z 2 ) 2 > q 2 /(si/z 2 ) and 
qf > q 2 > q 2 . The third parton is not generated, but its would-be-allowed space is marked by 
the leftmost almost-trapezoid line. 



Figure 6: The illustration of the kinematic of the BEV algorithm with n — 2 particles in the 
extreme case when the parton no. 1 is emitted with qf > so = s and the parton no. 2 is 
generated within the area not accessible for the parton no. 1 due to a higher IR boundary. 


space, is illustrated graphically in Fig. [6] It is, of course, an artefact of the BEV phase- 
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space parametrization, which in the FEV world corresponds to the phase-space reduction 
due to energy conservation. 

We omit here the discussion of the “kinematical cross-talk” between two parton show¬ 
ers, which means that for emissions with the common g 2 -ordering in two initial showers 
(as in any realistic PSMC), the emission in one shower reduces the available phase space 
in the other shower. This effect is easily incorporated in the kinematical construction 
(mappings) of PSMC. The only thing one has to watch out is the correctness of the soft- 
emission limit in the case of two and more emissions, see for instance the discussion in 
Ref. 


23 . This subject will be covered in a more detail in our future publications. 


3 KrkNLO methodology 


Very briefly, the essence of the KrkNLO prescription defined in Ref. 19 is that NLO 


corrections are introduced by a multiplicative weight on top of distributions from the LO 
PS acting, either inside PSMC or outside it, on a MC event record. This NLO weight is 
sensitive to the parton four-momentum with the highest q 2 


or maximum in the case 


of fey-ordering), although in Ref. 19 it was demonstrated that such a NLO multiplicative 


weight works also in the case of the angular ordering, provided that summation over all 
emitted partons is performed. 

In MCQNLO 1 21, the correcting weight is essentially additive and the NLO — LO cor¬ 
rection to PSMC is provided from outside in form of additional MC events, with a non¬ 
positive-definite weight. In parts of the phase space which are not covered by the LO 
PS, extra events provide the entire NLO distributions (positive weights), otherwise extra 
events, correcting the LO distributions to the NLO level, have typically (inconvenient) 
negative weights. 

In POWHEG H the entire NLO correction to LO PS is provided by an external 
MC generator - the highest fef, emission is isolated/subtracted from PSMC, following 
the double-logarithmic Sudakov exponential factor, and generated according to the NLO 
distribution outside PSMC, while trailing emissions with lower fe(suppressed by the 
Sudakov exponent) are left for generating within LO PSMC. 

Both KrkNLO and MCQNLO require good analytical control of the LO PS distributions, 
while for POWHEG it is less important. In addition, KrkNLO requires that LO PS fills in 
the entire NLO phase space with some non-zero distribution. The Sudakov suppression 
is also exploited in KrkNLO, but in a different way than in POWHEG. 

In the standard NLO corrections to a hard process with MS PDFs, part of the NLO 
corrections feature degenerate 1-dimensional longitudinal phase space with fe|, = 0 exactly. 
These corrections enforce in practice certain convolution on the MS PDFs, which in the 
case of POWHEG and MCONLO is done in-flight inside MC. In the KrkNLO prescription, 
the implementation of these corrections is moved outside MC (simplifying it). This leads 
to the use in KrkNLO of the modified LO PDFs, in the so-called MC factorization scheme. 
This reorganization is mandatory because one of the main aims of the KrkNLO method, 
that the NLO corrections are implemented with a well-behaved multiplicative positive 
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weight, is not compatible with such a degenerate collinear phase-space contribution to 
the NLO MC distributions. 

Last but not least, any scheme for correcting the hard process of LO PSMC to the 
NLO level requires a formal proof that the resulting distributions are indeed of the NLO 
class (no double-counting, no NLO leak). Such a proof in an algebraic form is not trivial, 
not only because the NLO total cross section has to be verified, but also any NLO- 
class observable (experimental event selections) has to be properly reproduced. In other 
words, it has to be done using functional space of all the NLO observables. In the case 
of KrkNLO such a proof was done in Ref. 19 , both algebraically and numerically, for 


a toy-model PSMC with the angular ordering. Here, we shall provide such an algebraic 
proof starting from the NLO-corrected multiparton distributions for realistic PSMC of 
the kind implemented in Sherpa and Herwig++ using the BEV algorithm. 

In the following, we are going to collect building blocks for the NLO weight, then we 
shall elaborate on the multiparton distributions of LO PSMC without and with the NLO 
weight of KrkNLO. We shall pay particular attention to the question of the completeness 
of the phase space in PSMC and to the equivalence between the backward and forward 
evolution algorithms in PSMC. Finally, we shall show that for an arbitrary NLO-class 
observable, KrkNLO gives the same result as simpler NLO calculations with the collinear 
PDFs, instead of PSMC, such as for instance MCFM 


29 


3.1 NLO-correcting weight 

Let us collect the ingredients for construction of the NLO corrections to the hard process 
of the Z-boson production and decay in the proton-proton collisions. 

The fully differential NLO cross section of the production and decay of the Z boson 
in the quark-antiquark annihilation process, with the simultaneous emission of a single 
real gluon, can be cast in a well-known compact form, see Ref. 


19 


d' , o"„ LO (a,P,Sl) = 


Cpa s dad/3 dip 


7r 


_ — dfi 
a(3 2i r 


da 0 {s, 9 F ) (1 - f3f dcr 0 (s, 0 B ) (1 - a) 5 


dfi 


di? 


, (3.1) 


where the Sudakov variable^ (a, f3) are spanned on momenta of the effective beams of 
q and q prior to the gluon emission, see Eq. (2.4). The Born differential cross section 

dao(s,d R ) £ ^ ,y „—, ,1 „.„ll ok j- 


exact expression) 


for Z-boson production and decay is well known (see for instance 


25 


The solid angle Q = ( 9 ,</>) is the direction of the lepton from the 


decaying Z boson in its rest frame and s = (1 — a — f3)si (see Fig. 


and 9 b depend on a and /? as well - their precise definition was given 6 in Ref. 30 


The angles 9 F 
The 


integration over luminosities of q and q is not yet included. 

How does the above compare with the distributions from PSMC also restricted to the 
single gluon emission? The differential cross section for the gluon emission from the quark 


5 For better readability of the formulae given in this section, we shall suppress “1” in the lower index 
of the Sudakov variables as well as other kinematical variables, such as s and qp B . 


In fact they coincide with the polar angles with respect to the effective beams defined in Eq. (2.12). 
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emitter (i.e. the ~ mp part in Eq. (2.6)) in PSMC reads 


,5 F , n n \ CfCXs dq2 F d V P / \ 7 Q\, n 

d a% q (a,P,i1) = ^Pqq{z)dz — {s,0) d f2 


C'fqi s dadfd dp l + (l — a — P) 2 dcr^ 
7T (a + /3)/3 27r 2 


(3.2) 


dfl 


(M), 


where qF is the evolution variable defined in Eq. (2.7) and 9 is another effective angle in 


Z decay specific to LO PSMC, for instance the so-called Collins Soper angle |31J. Adding 
the gluon emission from q simply amounts to the a <->• ft symmetrization, resulting in 


d 5 Vq°(a, P, M) = d 5 a% q + d b a B q = 

where 


Cpots dad(3 dp 1 + (1 — a — P ) 2 da 0 


7r a/3 27r 


—dfl 


dZl 


M), (3.3) 


j5 p _ ^ _LO _ j5 LO 

d — n d (J qq — vn F cl o = , 


qq 


ex. ~h /3 


d 5 a B = _—_ 

a ™ a + P 


LO 


d = m B d 


(3.4) 


The integration limits are not explicit, but they are the same as in Eq. (2.4). 


In our discussion, we shall often use the following objects: the additive NLO correction 


D, (2) = r" LO (a,>3, fi) - d?a£ (a, /?, fi) 


'5 „LO / 


and the NLO multiplicative weight for the qq channel 


wO( a p) = 1 + d5/3qq 

YV qq fj) L T 


<P°T 


(3.5) 


(3.6) 


The above weight is especially simple in the case of averaging over the angles in Z decay: 

2ap 


= i - 


1 + (1 — a — P) 2 1 


(3.7) 


and it can be used in approximate implementations of the NLO corrections. 
The analogous weight for the qg channel is 

(W (1) ) fl = = 1 + ^ + 2z "> . 

K q9 ,Q R (1 -z) 2 + z 2 


(3.8) 


After summing up the contributions from the two emitters, d 5 a is obtained, which 
is exactly the sameQ as in Ref. 19 . The important consequence of the above is that 
many details of the matching between the MS NLO corrections and the LO PSMC in 
the present KrkNLO implementation remains the same as in Ref. 19 . In particular, the 


In spite of the differences of the LO PSMC distributions for each emitter separately. 
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virtual plus soft-real correction, when the PDFs in the MC factorization scheme are used, 
is the same as in Ref. 19 and read^] 

-^mc = 1 + ^vsi ( 3 - 9 ) 

where 


Af s = -Ch 


2tt 


4 2 5 

-7T- 


A« s = 0 . 


(3.10) 


As one can see, this virtual+soft-real correction is constant, i.e. kinematics-independent. 


3.2 MC factorization scheme 

In this subsection we extend the definition of the MC factorization scheme, introduced in 


Ref. 19 for the quark-antiquark channel, to the NLO DY process with the quark-gluon 


initial state. For the completeness and convenience of the reader we first provide the main 
formulae for the qq initial state. 

The NLO gg-channel coefficient function for the DY process in the MS factorization 
scheme is given by fl9] 


crw = ?cvb(i-a(h 2 


2vr 


^ l + 


2 l±i! hzh! 

1 — z z 


(3.11) 


The corresponding coefficient function in the MC factorization scheme, defined in Ref. [19] , 
reads 


c? q c (z) = ^c F 


CXo 


S(l~z) ( ) — 2(! — z) 


Ota 


= 5(l-z)A vs --C F (l-z) 1 


(3-12) 


7r 


where Ays = Ay S is the virtual plus soft-real gluonstrahlung correction given in Eq. (3.10) 


From the above equations, following Ref. [19], we can obtain a qq contribution to the 
relation between the MC-scheme and MS-scheme quark (antiquark) PDFs: 


AC 2q (z) = - 


C% c { 


a s 


1±T hl (lvA + 1 


Similarly, for the NLO gg-channel contribution to the DY process we have: 


i^MS, 
g 1 

A) 


pMC, 

U 2 g 1 

A) 

a s m 1 . 

= — Tpt, — (1 

2vr 2 V 

AC 2 g( 

[Z) 

= C™(z) - 


(1 - z) s 


—z T 3 z T — 


iMC( 
'2 g \ 


a 


g T r ( [z 2 + (1 - z) 2 } In F_£lf + 2 2 (1 ~z)\. 


(3.13) 

(3.14) 

(3.15) 

(3.16) 


8 Note that B^ c , though conceptually similar, is not identical to the well known B used in the context 


of POWHEG 32 , as the quantity from Eq. (3.91 is defined in the MC scheme. 
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Using Eqs. (3.13) and (3.16), we can relate the MC-scheme quark (antiquark) PDFs 
to the corresponding MS-scheme PDFs in the following way 





AC2 q (z) + 



d z cMS 
.. J a 



AC 2l (z). 


(3.17) 

The above relation is universal, i.e. process-independent, at the NLO level. It is simply 
because it is defined uniquely with respect to the MS scheme. 

The gluon PDF is equal between the MC and MS schemes up to the 0(a s 2 ) corrections 
for processes with no gluons at the Born level, such as the Drell-Yan process considered 
in this work. Hence for the DY process we may use 


/rM )=/»(*.<?)• 


(3.18) 


As one can see, in the MC factorization scheme the NLO coefficient functions in both 
the qq and qg channels are substantially simpler than the corresponding ones in the MS 
scheme; in particular they are free of logarithmic singular terms. Since the latter terms 
are absorbed into the MC PDFs, i.e. they are included in a resummed way, one may also 
expect that the higher-order QCD corrections in the MC scheme are smaller than in the 
standard MS scheme. 

Let us finally add a comment on the universality (process-independence) of the MC 
factorization scheme. This issue was discussed quite extensively in Ref. [l9| for the DIS 
and DY processes, albeit restricting the problem to QED-like gluonstrahlung diagrams 
only. The main point is that the MC factorization scheme is in reality defined operationally 
as a modification of the MS scheme, consequently it inherits automatically the universality 
from the latteiJ^J On the other hand, it is true that the procedure of defining extra O(e 0 ) 
terms, added to the MS collincar counter-terms to define the PDFs in the MC scheme, is 
clearly guided by the inspection of a number of simplest physical processes, like DIS, DY 
and/or Higgs production. However, once this is done, these extra terms are frozen and 
the resulting (modified) counter-terms are ready to apply for any other process. They 
also define completely the PDFs in the MC scheme. What is new in the present work 
with respect to Ref. 19 is the inclusion of the quark-gluon transitions. Generally, the 


transformation of the PDFs from the MS to the MC scheme is a matrix in the flavor 
space. The DY process at the NLO level, discussed in the present work, Exes only a 
subset of terms in this matrix, while the remaining ones will be fixed by inspecting the 
NLO corrections to the process of gluon-gluon fusion into the Higgs bosoi: 9 10 At the next 


step, after including Higgs production in the game, the MC scheme will be fully defined 
and will be applicable to any process, including also more color particles in the final state. 


9 In particular, the MC scheme is neither the DIS nor the DY scheme, as can be seen from the fact 
that the coefficient functions for any of these processes in the MC scheme are not equal to <5(1 — z). 

10 They are set temporarily to zero in the present work, as for the DY process they become relevant 
only at NNLO. 
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3.3 Multiparton reality of PSMC 


In the above we have restricted ourselves to a single emission and hence omitted all the 
multi-emission reality of PSMC. Let us elaborate on that in a more detail now, because 
the NLO correcting weights are put not on top of single-emission distributions but on top 
of multiparton distributions of LO PSMC, so they have to be known and controlled for 
the hard process within the NLO precision, preferably in a closed algebraic form. 

Let us start from LO PSMC in the forward-evolution (FEV) formulation. The equiv¬ 
alent backward-evolution (BEV) formulation will be presented later on. Restricting our¬ 
selves temporarily to the pure gluonstrahlung case, the FEV differential cross sections of 


gluons emitted from the q and q emitters 11 * * reads as follows 


OO 


°mc = / dx F dx B df2 S / da: 


LO 


da. 


LO 


nn 

1=1 j =i 


d 3 P ?e, 


, 9 9 p~ S F^i-l^i) 

qf-i>g->qi c 



-S„(ql ,ql) -S„(ql ,ql) / n\ 1 n f ( „2 X F \ 1 ns ( „2 X B A 

F e s ~Tn^ sx F X B-> VV—-^mc Us’VV/V—^McOsi^F 

(IM Zj n F ' Zjnp ' Zj nB ' Zjub ' 


d Pj e 


-S B (qU’1i) 


(3.19) 


where Zi = YYi=i z h Z 0 = 1, q% = s and 7(z) = (1 - ^r)CP(^), 7{z) = P q q{z ). The 

principal evolution variable qf was introduced in Section |2| Eq. (2.7), and the labelling of 
the emissions starts from the hard process 12 , as in Fig. ST The emission distributions for 


the ladder labelled with p are the following: 


d 3 Pi = d*p? (s i:i ) = 


3 / 


dfdidati d(j>i - 


/ 3 i(cti + Pi) 2 vr 


(P(l Oli Pi) 0 ai> 0 d&i+l3i<l 


dqfdzi dpi 


e ( 


dqf dpi V(zi) 


Ti 


2 2 7T t7 ( 1_z <) 2s U > 9 t ? ^ > ( Z *) q2 2 p dZi \ 


Zi 


.^(1 -Zi) 2 Sij>q?J 


(3.20) 


where the Sudakov function reads 


S F (Qb,Qa) = S F( S ij\(lb,(la) = 


d Pi ( s ij)i 


(3.21) 


ql<q 2 i<qi 


and for the ladder b they look the same, except for the cq -H- Pi swap. 

The important variable s l0 entering d 3 pf ,B and S F,B for the single shower/ladder was 
already defined in Section [2] as s* = s/Zj, with s = sx F x B . For two showers, in any 
realistic PSMC, the emissions are generated (and the four-momenta are reconstructed) 

11 We adopt a convention in which d n = 1. 

12 This is unnatural for the present FEV scenario, but better suited for the BEV algorithm in the 

following. 
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simultaneously in both showers using the competition algorithm, in which a common q 2 - 
ordering in both showers is emerging in a natural watf^] Within such a common ordering 
for the two showers, the variable Sij = s/(Z[ Zf) includes all z 's from the emissions in 
both the ladders, starting from the hard procesepq 

Strictly speaking, the above implicit “kinematical coupling” of the two showers through 
variable prevents us from rewriting the distributions of Eq. (3.19), without any ap¬ 


proximation, into the traditional convolution of two LO PDFs and the hard process, as 
it was possible in the toy PS MC in Refs. (T9l 201. However, a slight modification of 


the kinematical coupling (modulo N 3 LO corrections) allows us to get from Eq. (3.19) the 
following equivalent factorized inclusive formula 


a 


mc = / dx F dx B dQ D^ c (s,x f )D 


da 


MC V'b %B ) 1^)1 


(3.22) 


where da/dQ is the hard cross section from Eq. (3.1) and the LO PDF D^ c is resulting 
from the FEV algorithm run separately for each single shower, written in form of the 
following time-ordered (T.O.) exponential 

00 p 

+ J2 (si) e - s ^ Sllq2 ^ ) d 3 p^{s 2 ) e~ s ^ q2 ' q2) ... 

n =1 J 

n 

x <*#•»**> n 

(3.23) 


i= 1 


S, X 


with Si = s/Zi, So — s, Qn = s. The other PDF, Dl 

F B /IK 

It is important to note that the objects D M ’ C appearing in Eq. (3.22) are not just scalar 


^MC 1 


,), is de 


ined analogously. 


functio ns but they have non-trivial and well-defined internal structure, as explicitly seen 
in Eq. (3.23). I 11 particular, the MC PDFs, , result from the Markovian process 
and thanks to kinematical mappings they respect the phase-space constraints exactly. 
Therefore, they are not equal to the standard DGLAP parton distributions functions, in 


particular they integrate emissions up the absolute phase-space limit, c.f. Eqs. (2.4) and 
(2.10), rather than stopping at some arbitrary scale Q 2 = / i 2 F . 


The reader may check, analyzing one and two emissions in a detail, that the effect of the 
above “kinematical coupling” of the two showers through the variable is conveniently 


absorbed by the construction of the four-momenta defined in Section 2.2 hence Eq. (3.22) 
is equivalent to Eq. (3.19) up to the N 2 LO level, i.e. neglecting the N 3 LO and higher 


corrections. The above equivalence could also be tested numerically, similarly as was done 
in Ref. fl9 . 


13 This method leads to forward-backward symmetric distributions, contrary to generating first the 
emissions from q and later on all the remaining emissions from q. 

14 In the PSMC jargon this is referred to as a recoil mechanism. 

15 See Ref. 33 for the precise formal derivation of the T.O. exponent from the Markovian FEV algo¬ 
rithm. 
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Why do we insist on the FEV representation of PSMC knowing that any typical PSMC 
is built using the BEV algorithm? The important reason is that in any methodology of 
combining fixed-order perturbative corrections with PSMC one has to make an algebraic 
contact with the standard diagrammatic perturbative calculations, including factoriza¬ 
tion theorems, resummation techniques, etc., which are all defined within the standard 
Lorentz-invariant phase space (LIPS). The FEV parton shower works directly within the 
LlPfj^J while in the BEV algorithm, the relation between LIPS and the PSMC distri¬ 
butions is obscured due to the presence of ratios of the PDFs in the BEV distributions. 
The way out is to define a framework in which the BEV and FEV distributions are by 
construction exactly identical. This requires certain non-trivial extra care and effort. We 
shall follow this path in the following. 

The PDF of Eq. (3.23) looks pretty standard, for instance, it is directly implementable 
in form of the FEV Markovian MC. Its peculiarity is already quite obvious for a single 
emission n — 1, where the upper phase-space boundary limit is not just regular qf < s, 
but more complicated q\ < (1 — Zi) 2 s/Z\, reflecting the realistic phase-space boundary of 
the hard process, while the lower boundary is just regular qf > q 2 , he. the same as from 
the ordering in PSMC. This peculiarity influences mainly the first emissions, closest to 
the hard processj^j as discussed in Section [ 2 } The consequence of the above peculiarity is 
that for our aim of the rigorous correspondence of the multiparton distributions between 
the FEV to BEV algorithms we need an auxiliary PDF with two competing factorization 
scales, Q 2 and s: 


D f mc (s\Q 2 ,x f ) = D F (s\q 2 ,x F )e s f( s \Q 2 ^) + ^ 


e ~S F (s\Q 2 ,qf) 


71=1 


x n {<?/>?(*) e-v<>*f.«f + i )e 02>3| _ i>Ji ,> j| j£i'(s| 9 , 2 ,x,) s xr=x . nj.,,, 


(3.24) 


2=1 


where q n+ 1 = q s and it coincides with the PDF of Eq. (3.19) for Q 2 = s = s/x F , he. 

d mc(^ x f) = d mc(s\s/x f , x f ) = L>m C (s|oo, x f ). (3.25) 


The main difference between this new PDF and the one introduced earlier in Eq. (3.23) 
is the 0Q 2 > q 2 1>q 2 >q 2 function, see Fig. 3j that restrict emissions only to the region below 
the scale Q 2 , which in principle is an arbitrary parameter, as long as Q 2 <C s. In that sense 


F.B 


the new PDFs, D^ c , are closer to the standard DGLAP parton distribution functions. 
With the above definitions we may work out the BEV algorithm providing the LO 


16 Modulo kinematical mappings, effective beam technique, recoils, etc. 

17 Already for the second emission, the limit < qf is more important than < (1 — Z 2 ) 2 s /( 21 ^ 2 ), 
which can be neglected modulo NNLO. 
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distributions, which by construction are exactly the same as from the FEV algorithm^} 

n F n B 

da ™ n B = ^ sxfXb n { d v n 

*=1 i=i ( 3 . 26 ) 

x e -A F ^ F \^ F ,^) e -A B ^ B \^ B ,^) DF ^ x ^ ZF(s,xJ dx F dx B di ?, 

ryF /'v* / V -f rpB /-v» / V B q Q r f' r P 

^ F / ^ j ^ B / ^ j 1 ^ odj F -ly B , 

where xf = UL=i z m includes z, from d 3 oof and similarly for x® from d 3 o;f. The single- 
emission distributions and form factors are defined as follows: 

dqfdzi dtfii 


d 3 u? = 


Qi 27T 


Kmc(^-i|^ 2 ) e-^cC^-d^d 

B^jiq{sx* f i x*/z) 


ifi-MC|A 17 ) — ^(1 — z) 2 sx* / z>q 2 -=r / , ^ . 

^Mc(z*|gLi,g, 2 ) = [ %- [ —K U c(x*\z,q 2 ), 


(3.27) 


n, 2 r A* 2 


and the shower label H — f, b was temporarily omitted wherever it was unambiguous to 
do so. 

Although it may not look obvious, the integrated cross section 


O T\ ,m — ^ ^ 

n F ,n B = 0 1 


L ° 

f MC 


da 


LO 


and all the multiparton distributions are exactly the same as in Eqs. (3.22) and (3.23). 
Let us present a sketchy proof of t half 3 } Directly from the BEV algorithm for the PDF 
of Eq. (3.23) of the showering quark we get 


)^ c (s,x F ) = D M c(sL|s,x)<J e 

I BEV 

ji 


-Awcis^Xp) 


F> + 


e f n f" 4 f - «-^^-)K l «(^,k 

^ L t= J<A Vi Jxi-X 


«?) 


D -A M c(xn\s,ql) 


(3.28) 


In order to prove that the above PDF is the same as from the FEV algorithm one exploit 


the fact that in the T.O. exponential representation of the double-scale PDF of Eq. (3.24) 


we may detach (factorize off) the first emission, obtaining the following integral equation: 
d mc4\ Q 2 i x ) = D F (s\q 2 s ,x)e- s r {slQ2 ’ q ° ) + 


r 4 /w 

Iq 2 T\ Jo Jo 


1 


dzi y{z 1 )9 [ i- zl) 2 S/zi>q 2 e SF(slQ2,q ^D MC (s/z 1 \q4x 1 )S x=ZlXl . 

(3.29) 


18 Modulo the kinematical coupling mentioned above. 

19 The result of the detailed proof is the same, of course. 
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Differentiating both sides by d/d{\wq 2 s ) and then integrating d(lnq 2 ), one obtains 


e -SMc(s\qtul) — c ~^mc (x\qh<ia) ^ MC 


DMc{s\ql,x) 


Duc(s\ql, x) 


(3.30) 


The above identity can be used many times in order to eliminate the ratios of the PDFs 
in the BEV formula and transform without any approximation the BEV distributions of 
Eq.(3.28) into the FEV distributions of Eq. (3.23). 


3.4 NLO weight in parton shower and its algebraic validation 

Having defined the framework in which the FEV and BEV distributions of LO PSMC 
are identical, we may define a MC weight correcting the hard process in the (notoriously) 
inefficient FEV LO PSMC, but having advantage of a straightforward connection to the 
perturbative expansion, and apply it within LO PSMC implemented using the efficient 
BEV algorithm. The NLO completeness of such a scheme can be proven analytically 
within FEV without any troubles, due to the absence of the ‘obscure’ ratios of the PDFs 
of BEV. The only extra cost is that in order to gain a perfect FEV -H- BEV compatibility, 
we had to introduce the special LO PDFs D(s\Q 2 ,x ) with two factorization scales, at 
least for the purpose of the discussion of the NLO completeness of the KrkNLO scheme, 
while in practice we shall be able to obtain D from the standard MS PDFs 

The differential cross section in which the hard process is NLO-corrected for the qq 
channel in the KrkNLO method reads as follow^} 

n F n B 

= (l + A vs + '£w$l(a*Jn + 'EwV(a‘J‘’))da™ nB , (3.31) 

i =1 j= 1 


where 


IE 11 - 1 = 
vv qq 


d% q d 5 


J qq 
d 5 (Tq? 


a; 


NLO 

qq 


d 5 ^ 


d 5 &q° 


= W$ - 1 , 


(3.32) 


see Eqs. (3.6) and (3.10). The structure of this weight, with the summation over the emit¬ 


ted gluons, is a straightforward generalization of the weight used in many MC programs 


with multiphoton exponentiated corrections, see for instance Ref. 34 
As pointed out in Ref. 


20 , in order to get the complete NLO corrections to the hard 


process, it is enough to retain in the sum in Eq. (3.31) only a single term, the one with the 


maximum k^. Since our g 2 -variable is practically identical to k^ we may retain only one 
term for a gluon with the maximum qf or qj , from one of the two showers. In the case 
of the BEV algorithm with competition, this gluon is just the one which was generated 
first. 

Similarly, as we have detached the first emission from the T.O. exponential represen¬ 


tation of the PDF in Eq. (3.29), we now factorize off explicitly the gluon emissions with 


20 


We adopt a convention that XX=i d n = 0. 
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the maximum q 2 and q 2 B from the product of two PDFs in Eq. (3.22): 


D mc( § ,x f )D^ c (s,x b ) =e S F^\ Sl F’is) e s b^b^s)D^ slq 2 , x F )Dl lc (s\q 2 s , x B ) 

+ [ d 3 pUs/zi F )e~ SF{slsi ^' q2l F ) e~ SBislsiB ’ q2l F ) —D^ c (s\ql F ,x F /z lF )D^ c (s\ql F ,x B ) 

J Zl F 

>1s«1i f < s i f 

+ [d 3 Pi {s/z lB )e~ Sp{s]si F' q2 B^e~ S B^ ba\ b )£> F MC ( s\q* , x F )—L>mc (s\q\ B ,x B /z± B ), 

J Z l B 

Qs«ll B <si B 


(3.33) 


where s i F = s/x F and si B = s/x B . 


In order to improve readability, we are going to omit temporarily (until Eq. (3.38)) the 

-p^F,B , . i 


argument s in the £)J^q(s| ...) and the S FjB (s\ ...) functions. Moreover, the “kinematical 
coupling” of the two showers is will be restored, that is we are going to replace in the 
following Si F , Si F —> s = s/(x F x B ). The above decomposition is easily exploited in 


calculating the NLO cross section with the MC weight of Eq. (3.31) truncated to a single 
term with the maximum q 2 : 


a 


NLO 


OO n n 

E /<?.=/ 


da ntn„ = I dx F dx B dfi ^{sx F x B ,0){l +A vs )e s F( s ^) e 


^ -^MC ids i j *^b ) d LO (x F , X B ) 


dx F dx B dfi 


d 3 pf(sx F x B )(l + A vs + WA{ki))dd-(sx F x B zi F ,6) 


l 1 ] i 
qq 


da 


dfi 


'q 2 s<ql F <s 


x e -S F {s,qi F ) e -S B {s,ql F ) n F f 


F -^Mc( ( ?li?> {dl F i X B ) AmLO F > X B , Zip, kyi *) 


+ 


d 3 pf(sx F x B )(l + A vs + (ki))^-(sx F x B zi B , 9) 


qq 


da 


dfl 


1s<1l B < S 


xe s f^I b ) 6 S B^l B ^Dli C {ql x F )D^ c {ql x B )J n vo{x F ,x B ,Zi B ,kl T ) >, 


(3.34) 


where we have introduced x F and x B before the emission, as in a typical fixed-order NLO 


calculation. They are related to xp and xb from Eq. (3.33) in such a way that if the 


emission happens on the leg F, then xf = x F /zi F and x B — x B . And symmetrically if 
emission occurs on the leg B. The same expression we can get directly from distributions 
of the BEV algorithm, as shown explicitly in Appendix [Aj 

For the purpose of the proof of compatibility of the above formula with the fixed- 
order calculation (MCFM) for the entire class of the LO and NLO observables, we have 
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introduced in the above the jet functions , Jlo(x f ,x b ) and Jnlo( x f , x B i z i> ^it)j which 
satisfy the following properties: 

JNLo( x F ,x B ,zi,kl T ) ->• Jlo(x f ,x b ) for k\ T -> 0, or zi ->• 1, (3.35) 


often referred to as an infra-red safety requirement. In the experimental practice, the 
function Jnlo represents the most general 2-dimensional histogramming in the transverse 
momentum of the gluon (or the heavy boson) and the z\ variable, related to the rapidity 
difference between the heavy boson and the emitted gluon (or the ratio of the effective 
mass of the heavy boson and the heavy boson plus gluon system). The meaning of the 
above limit k\ T —^ 0 is that for the first bin in k 2 T , the one including the k\ T = 0 point, 
we are not allowed to do any binning in z\ (we have to sum up inclusively over all z i, 
similarly as for Jlo we sum up inclusively over all k\ T ). 

In order to recover the fixed-order NLO formula of MCFM, all terms 0(a 2 ) have to 
be carefully eliminated. It is easy to do it for the ~ d 3 p Wj^ parts, which are formally 
0(a\). We replace in them -Dm C (s|<2i f i s , X F ) —>■ -Dmc(^I s > x f ) and also eliminate safely 
the Sudakov exponent^ 1 } 


cr NLO [J] = J dx F dx B dn |e S BMD^ c ( q 2 s ,x F )D^ c ( q 2 ,x B ) 

da 


x (1 + A vs )—(s 1 ,6)J LO (x F ,x B ) 


dx F dx B dft 




d Pi (si)e 


^ BKa,HlF> D MC ( qi x F )D MC ( qi x B ) 




+ d 3 pf( Sl )e S F {s ^ B ) e Ss{s,q2lB) Dli C { q l x F )Dli C { q l x B ) 


<ls«ll B <S 


da 


x (1 + A vs )—(sx F x B zi,d)J NLO (x F ,x B ,zi,k 1T ) 


da 


+ dx F dx B di1 / (d 3 pf(s 1 ) + d 3 pf(s 1 ))WA(k 1 ) — (sx F x B Zi,9) 


o«?i F ,gf s <s 


^ X F )"^MC > X B ) '^NLO ( X F > x b 5 Z\i kuj , 


(3.36) 


where we were also allowed to replace zi F ,z lg —* z i, because in terms of the Sudakov 
variables before the emission it is the same variable. We have also introduced si = sx F x B . 

The remaining parts ~ (1 + Ays) are less trivial. In order to use again the identity 
of Eq. (3.33) for folding in three integrals within the {... } braces back into a simple 


21 The collinear singularity in q 2 —» 0 is killed by . so f dq 2 /q 2 e 
contribution. 


cannot give a ~ 1/c 
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product of two PDFs, we have to do something with the Zi-dependence in da/dft and 
the fc^-dependence of the jet function Jnlo in the second part. The solution is to add 
and subtract a term proportional j^(si,9)J LO (x F ,x B ) in the second part, regroup and 


use again Eq. (3.33): 


/ d(j 

dx F dx B dfl PM C ( s - x f) Z} Mc( s . ;i: fl)( 1 + A V s)-^(s 1 ,6)J lo (x f ,x b ) 
dx F dx B di / d 3 p^( Sl )D^ c (ql F ,x F )D^ c (q 2 x B ) 


+ / d 3 pf^D^ql^x^D^ql^x^ 


(3.37) 


x 


^( SXFX B z lJ) J NL°( x F , x B , z l,kl T ) - ^(s u 6)Jtx)(x f ,x b ) 


dx F dx B dfl / d 5 p^9 fE'jy(/ci)DM C (s,a: F )DM C (s,a: s )J NL o(x F ,a: s ,^i,/ci T ), 


where we have profited from finiteness of the integrals, in order to omit (1 + Ays) an d 
the Sudakov exponents wherever possible and to set the lower integration limits of q 2 to 
zero. 

The final clean-up of 0(a 2 ) terms involves the change q\ —» s and q\ —> s in the 
PDFs in the second integral and recombining it with the third one: 


a 


NLO 


/I 

dx F dx B dfl (1 + A vs ) — (s 1 ,e)J LO (x F ,x B )D^ /lc (s\s,x F )D^ AC (s\s,x 1: 
dx F dx B di1 / \d 5 p^° [1 + W^(k 1 )]J NLO (x F ,x B ,z 1 ,kl T ) 


dcr 


— (d Pi + d pf)—( Sl ,6)J LO (x F: x B ) fD MC (s|s, x f )D mc (s|s, x b ). 


(3.38) 


We may finally eliminate the MC weight of Eq. (3.32) going back to the NLO distributions: 


da 


<t nlo [J] = / dx F dx B di1 (1 +A vs )^(s ll 9)J LO (x F ,x B )D^ c (s,x F )D^ c (s,x B ) 


+ I dx F dx B dfi\d 5 pgg Jnlo (x f ,x b , z h k( T )-d 5 p^ Jlo(x f ,x b ) 

X ^Mc(®> S f)-^Mc(^ X b)i 


(3.39) 


where (s, x) of Eq. (3.25) was also recovered. 

The above looks like an example of the fixed-order NLO calculation formula employing 
the technique of soft-collinear counter-terms following the Catani-Seymour (CS) work 1241, 
with the explicit definition of an arbitrary NLO observable using the J-function. For the 
sake of completeness, such a formula in the standard MS scheme using the CS method is 
shown explicitly in Eq. (B.7) in Appendix [B] for the same qq channel. 
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However, t here a re two important differences between the formulas in Eq. ( B.7 ) and in 
the above Eq. (3.39): missing the ~ 5(kf T )£ term and non-MS PDFs in Eq. (3.39). The 
real emission integral with subtraction, in spite of slightly different notation, is identical. 
These differences are, of course, due to the differences between the MC and MS factoriza¬ 
tion schemes, and are well understood, see the discussion in the following Section [4] and 
in Appendix [B] The important bonus of the comparison of Eqs. (B.7) and ( 3.39[ ) is that 
we can establish the relation between the quark PDF used in the parton shower MC and 
the corresponding PDF of the MS scheme, modulo NNLO, in a solid and unambiguous 
way: 


f„ (Q ,x) = D mc (Q ,x) = D mc (Q I Q /x,x), 


(3.40) 


where is the quark PDF in the MC scheme defined in Eq. (3.17). 

In this way we have proven algebraically that the KrkNLO scheme is equivalent to the 
fixed-order NLO calculation in the entire functional space of the NLO-class experimental 
observable^] 


3.5 Summarizing the KrkNLO method 

In the following, we summarize and comment on the key elements of the KrkNLO method, 
starting from the qq channel only and then adding the qg channel. 


1. The first element is the reorganization of the NLO corrections, the same as Ref. 
Starting from the standard MS calculation, a strictly collincar part of the NLO 
corrections is removed thanks to the cle ver re definition of the LO PDFs from the 
MS scheme to the MC scheme, see Eq. (3.17). In a nutshell, in the standard MS 
procedure, one subtracts from the diagrammatic real + virtual results a pure ~ 
P(z)/e pole, while in the case of KrkNLO one subtracts the r M c(z,e) function, 
which contains extra non-pole terms. As discussed in Ref. 19 , see also Eq. (B.ll), 


this r MC function is just the integral over the LO distribution of the single-gluon 
emission in the parton shower MCp^j written in 4 + e dimension, together with the 
Sudakov form factor, obeying typical parton shower sum rules. This element of the 
KrkNLO prescription remains the same as in Ref. 19 , provided the LO distribution 
is the same as in the present study. This assumption is valid because the LO 
distribution of Eq. (3.3) is indeed the same in Ref. 19 and in the Catani-Seymour 
inspired implementations of the single-gluon emissions in Sherpa and Herwig++. 
What remains to be checked is whether the upper phase-space limit, a + fi < 1, is 
not spoiled in a given PSMC by the use of the backward evolution algorithm. We 
are going to come back to this point shortly. 


2. The second clement of the KrkNLO prescription is the construction, implementa¬ 
tion, and algebraic validation of the multiplicative weight introducing the NLO 

22 It seems that this kind of an explicit rigorous, albeit tedious, algebraic proof is not available for the 
other methods of combining the NLO corrections with the LO parton shower. 

23 It is the so-called soft-collinear counter-term. Although it is defined with the help of PSMC for the 
DY process, its definition is universal. 
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corrections in the multiparton environment of PSMC. The weight of Eq. (3.6) is 


implemented following Eq. (3.31). The summation over gluons is necessary only for 
the angular ordering, while for the g 2 -ordering or fey-ordering it is enough to keep 
just one term for the gluon next to the hard process. Validity of above method can 
be proven rigorously, see Section |3.4[ in the case when the initial-state PSMC is 


realized using the forward Markovian evolution (FEV) algorithm, in the presence 


of any observable-defining function J, see Eq. (3.35). Such a FEV algorithm would 


be terribly inefficient for the resonant Z-boson production, but it is mathemati¬ 
cally perfectly well defined. (In such a gedanken FEV scenario with the complete 
phase-space coverage, the output PDF from the MC at the hard process scale will 


be automatically in the MC scheme, following Eq. (3.17), provided that the input 
PDF at low q 2 s is in the MS schemcF 4 !) 


3. The above validation proof of the NLO weight has to be extended to the (efficient) 
backward evolution (BEV) scenario of the typical PSMC, such as Sherpa or Her- 
wigTT- It was shown in Section 3.3 that this is feasible, by means of formulating 


the twin FEV and BEV algorithms, using the same PDFs, which produce exactly 
the same multiparton exclusive distributions. What is highly non-trivial is that the 
full phase-space coverage is not lost. The p rice for this was that we had to introduce 
the auxiliary PDF D(s\Q 2 ,x) of Eq. (3.24) with two competing factorization scales, 


in addition to the single-scale PDF of Eq. (3.23). Luckily, at the end this auxiliary 


PDF can be finally eliminated from the BEV algorithm. Its main role is to provide 
theoretical control for the FEV —> BEV transition, desired in view of the previous 
point. 


4. The elimination of the auxiliary parton distribution function D(s\Q 2 , x) can be seen 
best with the analysis presented in Appendix [A] where the transition FEV —>- BEV 
is analyzed in a fine detail for n — 0,1 emissions. In E q. (A.l) for n = 0, we see that 
we may replace D(s|s,a:) —> D(s,x ) using Eq. (3.25). For n = 1 in Eq. (A.6), one 
may show that, modulo (9(a 2 ) terms, the ratio of D-functions can be replaced with 
any kind of the LO PDFs and the exponents, exp(— A), can be neglected. What 
cannot be approximated there, it is the exact implementation of the evolution kernel, 
keeping the correct phase-space limits due to the ^-functions. This is, however, not 
the problem, because any standard BEV implementation with the veto algorithm 

We conclude that, in the BEV scenario, all D-functions 


does that correctly 
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go away and only the single-scale PDFs in the MC factorization scheme seen in 
Eq. (A.l) are left finally in the game 2 ®' 


24 In other words, such a FEV MC perforins not only the LO DGLAP evolution but also the transition 
from one to another factorization scheme. 

25 More detailed analysis of this kind was done but is left beyond the scope of the present paper. 

26 One can show algebraically that the effective PDF generated by the above BEV at low g 2 will be in 
the MS scheme - the BEV algorithm is not only undoing the LO evolution, but also the transition from 
the MS to MC scheme. 
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MCFM 

Sherpa 

Herwig±± 

* 3 tot [pb] 

936.9 ±0.1 

937.2 ±0.2 

937.0 ±0.6 


Table 1: Values of the total cross section with statistical errors at the Born level for the 
Drell-Yan process in the MS scheme. 


In the above analysis we have omitted qg channel. For the Z/ 7 *-boson production 
process however, is not difficult to include it following the same steps described above} 27 ] 
In our implementation of the KrkNLO method in PSMC, whenever in the first step of the 
BEV PS algorithm the transition from the quark (antiquark) to the gluon is generated, 
we associate with such an event the weight computed according to Eqs. (3.8) and (3.9). 
This weight corresponds to observables averaged over the Z/ 7 *-boson decay angles and 
is sufficient for the current paper, as we are going to concentrate on the Z/ 7 *-boson 
transverse momentum and rapidity distributions. The discussion of the Z-decay leptonic 
observables, as well as more details on the implementation of the qg channel, is reserved 
for a separate publication. 


4 Fixed-order NLO benchmarks 

For the numerical evaluation of the cross section^] at the LHC for the proton-proton 
collision energy of yfs = 8 TeV we have chosen the following set of Standard Model (SM) 
input parameters: 


M z = 91.1876 GeV, r z = 2.4952 GeV, 

M w = 80.4030 GeV, r w = 2.1240 GeV, (4.1) 

G „ = 1.16637 x 10 “ 5 GeV“ 2 , m t = 173.2 GeV, 
a s (M 2 z ) = 0.13938690, 

and the G^-scheme |35] for the electroweak sector of the Standard Model. To compute 
the hadronic cross section we also use the MSTW2008 LO set of parton distribution func¬ 
tions |36], and take the renormalization and factorization scales to be /jl 2 r — /jl 2 f — M z . 
The only detector acceptance cuts are imposed on the invariant mass of the final-state 
lepton pair (Z/ 7 *-boson): 

50 GeV < Mu < 150 GeV. (4.2) 

In order to check that our settings are identical in all used programs, we began with 
the comparisons at the Born level. The results presented in Table, [l] show a very good 
agreement (within statistical errors) between different programs. 

2 'All ingredient distributions are well known in the literature. 

28 Unless stated otherwise in the text. 
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Figure 7: Ratios of MC PDFs to the standard MS PDFs at Q 2 = 100 GeV. The two curves 
in each plot correspond to the case with (red-solid line) and without (dotted-blue line) the 
convolution with the gluon distribution in Eq (3.17). 


4.1 PDFs in MC scheme 


In Section 3.2, we have introduced the MC factorization scheme and explained why this 


scheme is better suited for matching the NLO results with the parton shower than the 
standard MS scheme. 

The MC factorization scheme comes with a new set of parton distribution fun ctions 
that can be obtained from the standard MS PDFs via the relations given in Eqs. (3.17) 


and (3.18). Note that the convolution terms in Eq. (3.17) introduce the dependence on 


the renormalization scale via the strong coupling multiplying the expressions for AC^q 
and AC 2 g in Eqs. (3.13) and (3.16). In what follows, we shall always choose that scale to 
be equal to the scale Q 2 = fi 2 F of the MS PDF. 

For clarity of the discussion, before presenting the complete NLO corrections to the 
Drcll-Yan process, we shall be starting in what follows from results limited to the pure 


qq channel. This case corresponds to setting the gluon coefficient functions of Eqs. (3.14) 


and (3.15) to zero. That in turn implies that, for the two schemes to be consistent at 


0(a s ), the last term in the relation between the quark PDFs in the MS and MC schemes 


from Eq. (3.17) should be omitted. 


set 


The MC PDFs for light quarks and light antiquarks, obtained from the MSTW2008 LO 
shown in Fig. [7] through ratios to the MS PDFs. The curves correspond to 


36 


are 
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the choice of Q 2 = 100 GeV. The same scale was used as the argument of a s . Each panel 
shown in Fig. [7] contains two curves. The solid line anticipates including both the qq and 
qg channels, hence it was obtained exactly with the formula (3.17). On the other hand, 
the dotted line corresponds to the case of the pure qq channel, hence only the convolution 
with AC 2q was used, as explained earlier. 

As see n in Fig. [7J the MC PDFs without the gluon, e.g. without the last term in 
Eq. (3.17), are very similar to the MS quarks at low and moderate x. The MC PDFs in 
this region are only about 2% higher. At large x the ratio increases, reaching the value 
2-3 near x — 1. Note that the MC PDFs are always higher than the standard MS quark 
PDFs, which results from the fact that the second term on the RHS of Eq. (3.17) is always 
positive. Adding the gluon component to the MC quark makes a little difference at large 
x, whereas the small and moderate x regions are affected significantly. We see that the 
convolution with AC 2g is negative and, at small x, it can reduce the MS quark PDFs even 
by 20%. We have found a very similar picture for other quark flavors as well for other Q 2 
values. 

Before moving to NLO, it is interesting to check how the LO result changes when 
switching from the MS to MC PDFs. From the LO point of view, both PDFs are equiva¬ 
lent, as the differences are formally of the 0(a s 2 ) order and higher. We can see, however, 
some numerical differences in the LO distributions between those cases. And, indeed, 
as shown in Fig. [8j the differential distributions of the dilepton mass and the dilepton 
rapidity differ at LO depending on whether the MS or the MC PDFs are used. In the 
case of the MC PDFs without the gluon convolution, c.f. Eq. (3.17), the mass spectrum 
is ~ 5% higher than that with the MS PDFs over a broad range around the peak. On 


the other hand, when we include the gluon convolution, the result with MC PDFs gets 
up to ~ 20% below that of MS. The dilepton rapidity distribution with the MC PDFs 
with (without) the gluon is smaller (larger) by a similar amount for central rapidities and 
grows above 100% in t 
x behavior seen in Fig. 7 

Here, X\ and x 2 are the usual hadron’s energy fractions of the incoming partons, while 
yz is the rapidity of the produced Z-boson. The large rapidities in the forward or back¬ 
ward direction correspond to one of the x’es being large and that is the region where the 
differences between the MS PDFs and the MC PDFs are substantial. 


re forward and backward regions. This can be related to the large 

via the formula X\ 2 = exp (Ayz) valid for 2 —>• 1 kinematics. 

Vs 


4.2 MC and MS schemes at NLO 


The cross section for the Drcll-Yan process at NLO in the MS scheme can be schematically 
written as 


^dy = 



z) + ot s C q 


MS 


®fj(x 2 ) 


^fi(x 1 ) (8) a s C™ s ® fj(x 2 ), (4.3) 


h3=Q,g 
*7V' 


i=q,q;j=9_ 
1 = 9 \j=q,q 


where fi(xi t2 ) are the standard MS parton distribution functions, whose scale dependence 
is understood, and <E> denotes the convolution via integration over the z variable. The 













m z [GeV] y z 


Figure 8 : Distributions of the transverse momentum and rapidity from the LO calculation 
from the MCFM program with the standard MS MSTW2008LO PDFS and with PDFs in the 


MC factorization scale with and without including the gluon convolution in Eq (3.17). 


first term in the square brackets is just the Born contribution, which is followed by the 
NLO correction in the qq channel. The second sum corresponds to the NLO contribution 
from the qg channel. In both cases, the functions C ) 333 have non-trivial ^-dependence 
(suppressed here for cl arity of notation) and they are related to the coefficient functions 
from Eqs. (3.11) and (3.14) as Ci)l s = A- CM 1 | (z). Hence, the dependence on a s is 


explicit in Eq. (4.3) and all the following formulae of this section. 


When going to the MC scheme, the MS coefficient functions in Eq. (4.3) need to be 


replaced by their MC scheme counterparts and the PDFs need to be transferred to the 


MC scheme as well according to Eq. (3.17). This leads to the formula 


a 


MC 

DY 


^ (fi(x i) + fi(x i) (8) a s AC q + f g (x i) (8) a s AC i 




5(1 — z) + a s C ( 


MC 

q 


(x 2 ) + fj(x 2 ) <8> a s AC q + f g (x 2 ) 18 ) a s AC & 


+ y (fi(x 1 ) + f%{x 1 ) <8) a s AC q + f g (x 1 ) 18) a s ACg) ® a s Cf c 0 f g {x 2 ) 


i=q,q 


fg(x 1 ) 0 a s cy 0 ( fj(x 2 ) + fj(x 2 ) 0 u s AC q + f g {x 2 ) 0 a s AC g ) , (4.4) 


j=q,q 


where, again C^^ g (z), with the MC coefficient functions defined in Eqs. (3.12) 


and (3.16). 


and (3.15), and similarly for AC qjg = A- AC 2qt 2 g(z), with the latter given in Eqs. (3.13) 
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LO 

00 y [pb] qq channel 

°dy [pb] both channels 

MS 

936.79 ± 0.30 

MS 

988.9 ±0.30 

MC 

989.18 ±0.32 

MC 

778.8 ± 0.20 

in which: 


in which: 


0(a s ) 

fqfq 

25.79 ± 0.04 

fqfq ± fqfq 

54.8 ± 0.4 

fqfq 

25.79 ± 0.02 

fqfq ± fqfg 

-271.4 ±0.4 

o 

(« 2 ) 


0.64 ±0.01 


6.70 ±0.20 


Table 2: Values of the total cross section with statistical errors for the the Drell-Yan process 
at LO in MS and MC factorization schemes. The results were obtained with MCFM 6.6 1291. 
The 0(a s ) admixture in the MC results is split into contributions proportional to various terms 
of Eq. (4.4). 


By construction ( c.f . section 3.2), terms proportional to a given partonic luminosity, 
fifj, are identical in both factorization schemes up to the order 0(a s ). The results in 
Eqs. (4.3) and (4.4) differ however at 0(a s 2 ), which is beyond NLO and therefore such a 
difference is allowed. 

In order to validate our transformation from the MS to the MC scheme, we have 
performed an explicit check of the equivalence of the two schemes up to 0(a s ), as well 
studied numerical importance of the higher order terms. The calculations were performed 
using our standard setup defined at the beginning of Section [4j 

Let us start from comparing the cr oss s ection s at LO. This corresponds to setting the 
coefficient functions C (hfJ = 0 in Eqs. (4.3) and (4.4). As shown in the previous section, 
the LO cross sections will differ between the two schemes because of the PDFs. To check 
the extent to which this happens, we performed an explicit computation with MCFM |29j 
in both factorization schemes for either the pure qq channel or both channels. 

Indeed, as demonstrated in Table[2j the LO MS and MC cross sections are not identical, 
both in the case of the pure qq channel, where the LO result with MC PDFs is ~5% higher, 
and in the case where both channels are included, where the MC result is ~20% lower. 
The beyond-LO, 0(a s ) and 0(a s 2 ), terms are also given in the table. We see that most 
of the difference comes from 0(a s ). The 0(a s 2 ) terms are in fact very small, below 
1% in both cases. It is interesting to note that the difference between a^y bi the two 
schemes, for the case with both channels, comes primarily from the term proportional to 
a s f q ® AC g 0 f g , hence it originates from the large gluon luminosity. 

We turn now to a similar comparison at NLO. Table [3] gives the NLO-only results 
in the two schemes for the cases of qq and both channels. At NLO, we need to be 
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NLO 

°dy [pb] QQ channel 

ctqy [pb] both channels 

MS 

336.36 ±0.09 

MS 

157.9 ±0.10 

MC 

352.96 ±0.09 

MC 

305.8 ±0.10 

in which: 


in which: 


0(a s ) 

fqACqfq 

25.79 ± 0.04 

f q AC q fq ± fqACfq 

355.00 ±0.29 

fqACqfq 

25.79 ±0.02 

f /iMC f i f f 

Jq^q Jq ' Jq^q Jq 

-17.10 ±0.14 

f (-iMC f 
Jq,q^q J q,q 

284.77 ±0.08 

f Omc f 1 f A mc f 
Jq^g h ' Jq^g h 

-180.10 ±0.10 

sum 

336.35 ±0.09 

sum 

157.80 ±0.34 

0(a£) ± 0(cts) 

MC 

16.61 ±0.02 

MC 

147.9 ± 0.20 


Table 3: Values of the NLO contribution to the total cross section with statistical errors 
for the the Drell-Yan process in the MS and MC factorization schemes. The results were 
obtained with MCFM 6.6 1291 and its version adjusted to the MC scheme, MCFM*. The 
0(a s ) corrections in the MC results is split into contributions proportional to various terms of 
Eq. gg. 
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careful what we take for the coefficient functions. In the MS scheme, we just use the 
original MCFM 6.6 29 implementation. However, in order to carry out consistent NLO 
calculations in the MC scheme, we had to change the coefficient functions from C^ s 2g (z) 
to C%%{z). We dubbed this modified version of the program MCFM* and used it together 
with the MC PDFs to compute our NLO predictions in that scheme. 

We see that also the subleading corrections to the Drell-Yan process are not the same 
in the MC and MS schemes (numbers in the two first lines of Table [ 3 J. This is allowed 
provided that all the difference comes from terms of the order a s 2 and beyond. As shown 
in Table [3j by extracting only the 0(a s ) terms of the MC result and summing them 
up, we recover the MS cross section exactly. There, we also give subleading correction 
introduced via the MC PDFs and we see that they contribute at most 5% to the NLO 
correction in the MC scheme in the case of the pure qq channel but can be quite sizable 
when both channels are considered. 

Finally, it would be interesting to fit the quark and gluon PDFs in the MC factorization 
scheme directly to experimental data, in order to minimize higher-order effects. For that 
we would need to define the gluon PDF in the MC factorization scheme using a process 
in which the gluon PDF enters already at the LO level, e.g. the Higgs-boson production 
process. This will be done in our future publication. 

Let us conclude that we have successfully validated the MC factorization scheme by 
explicitly showing that all numerical differences w.r.t. the MS results are of the order of 
a s 2 and a s 3 . We emphasize that the validation of the new MC scheme presented in this 
section provides a highly non-trivial check as various components, given schematically in 


Eq. (3.17), come from different parts of the calculation (PDFs, coefficient functions) and 
the agreement up to a s is much more sophisticated than a purely algebraic relation. 


5 Results for NLO with parton shower 

Numerical implementation of KrkNLO is presently done mainly using the Sherpa PSMC 

^with the dipole organization of the parton shower distributions 
see Ref. 


of Ref. 11 version 2.0. 


inspired by the Catani-Seymour work 24 , see Ref. 37 . However, the evolution variable 
in Ref. |37| is chosen to be the transverse momentum distribution, while in the actual 


Sherpa 2.0.0 implementation it is the q 2 ~ a(a + (3) variable of Section [ 2 J see Ref. 38 . 
The dipole shower implemented in F!erwig++, see Ref. 1231, is quite similar. Actually, the 
choice of the evolution variable is not critical in the KrkNLO method, as long as the full 
coverage of the phase space for the emission closest to the hard process is assured, and 
the PSMC distribution is under the perfect control. In fact, for any choices of the dipole 
PS evolution variable in Sherpa or Fierwig++ the same distribution of the single-gluon 
emission is obtained, provided that the contributions for the quark and antiquark emitters 
are added. For the Sherpa MC program, we have checked numerically the completeness 
of the coverage of the phase space of the first emission from the backward evolution by 
examining the gluon distribution on the a, j3 plane. 
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https://sherpa.hepforge.org/doc/SHERPA-MC-2.0.0.html 
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In the following, the KrkNLO-matched results will be compared mostly with the results 
of the MCONLO technique, implemented on top of the same parton shower within SherpaJ^J 
Using the same parton shower in both cases ensures that all differences in the matched 
results come purely from the differences in the methods themselves rather than details 
of the shower implementation. Nevertheless, the difference between the dipole showers 
currently implemented in Sherpa and Herwig++ turns out to be small. From the point 
of view of our study, it amounts mainly to using the q variable in the former and kx 
in the latter. These two evolution variables are, however, very close to each other in 
practice. Hence, in our final discussion, we shall also compare the KrkNLO results with 
those of POWHEG, as implemented in Herwig+Up 1 ] using an automated setup based on 
the Matchbox framework 


39 with adaptive sampling |40|. 


In the numerical validation of the KrkNLO approach, we shall first compare it to 
the fixed-order NLO calculation, as implemented in the MCFM program. Later on, the 
comparison with the MCQNLO implementation of Sherpa will be discussed. In order to 
minimize the influence of various higher-order effects on the difference between KrkNLO 
and MCFM, we shall initially limit our numerical exercises to a simplified version, in 
which only the qq channel is kept. Another initial limitation, with the same aim in mind, 
will be to profit from the option in Sherpa to stop the backward evolution in the parton 
shower just after the first emission (starting from the hard process). This will help us to 
eliminate the influence of the differences in the so-called recoil schemes 41 , which start 


to play a role from the second emission onwards. The running ay will be also introduced 
gradually. Since we study the differences in the matching methods, we have switched off 
in Sherpa non-perturbative effects including the intrinsic kx, multiparton interactions and 
hadronization. However, we would like to stress that there is no problem with switching 
them on for the KrkNLO method. 

In all the following results for the KrkNLO method the PDFs in the MC factorization 


scheme, discussed in Sec. |4.1[ will be used. The factorization scale will be set to /ip = M z . 
As for the choice of the renormalization scale in the matched NLO + PS results, we note 
that there is always some level of arbitrariness. In the pure fixed-order NLO calculation, 
in the case of DY, one usually sets /j,r = /ip = M z . On the other hand, the LO parton 
shower uses p 2 R = q 2 , with the latter being closely related to the evolution variable of the 
shower. The parton shower is unitary, hence the choice of the argument of ay (/ip) does 
not influence the total cross section, it changes, however, shapes of some distributions, 
e.g. of px.z- The differences between a s (q 2 ) and a s (Mf) are formally of the higher order, 
that is beyond NLO and, from that perspective, they are equivalent in the context of 
NLO + PS matching. They can, however, be numerically relevant. 

In what follows, for the KrkNLO method we adopt the following procedure. The 
argument of ay in the virtual correction is always set to Mf. For the real correction, all 
emissions except for the first one (in BEV) are set by the parton shower according to its 
current q 2 value. For the first emissions, we consider two choices: q 2 and Mf, both in the 


30 Unfortunately, the POWHEG method is not available in Sherpa 2.0.0 framework. 

31 For our comparisons we used the version of POWHEG without restrictions of the phase space, which 


is similar to the original POWHEG method described in Ref. [3 . 
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°to t [pb] 

MCFM 

1273.4 ±0.1 

MCONLO 

1273.4 ±0.1 

POWHEG 

1272.1 ±0.7 

KrkNLO a s (q 2 ) 

1282.6 ±0.2 

KrkNLO a s (M 2 z ) 

1285.3 ±0.2 


Table 4: Values of the total cross section with statistical errors for the Drell-Yan process (the 
qq channel only) from the KrkNLO method compared to the fixed-order result of MCFM and 
the results of MCONLO and POWHEG. 

8 TeV: qq channel (1st emission only) 8 TeV: qq channel (1st emission only) 



Pt,z [GeV] y z 

Figure 9: Comparisons of the transverse momentum and rapidity distributions from MCFM, 
MCQNLO and two versions of KrkNLO, one with fixed a s (M 2 ) (also in MCFM) and another 
one with running a s depending on the transverse momentum squared. KrkNLO method is 
implemented within Sherpa parton shower MC. The exercise is restricted to the qq channel 
only and the parton shower backward evolution is stopped after the first emission, the closest 
to the hard process. 


real correction and in the Sudakov form factor, to keep the shower unitary. The difference 
between the results corresponding to these two choices will be indicative of the size of 
beyond-NLO terms. 
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5.1 Initial results for qq channel only 

The first round of numerical tests of KrkNLO focuses on comparisons with MCFM and 
MCONLO for the qq channel only and stopping the parton shower after the first emission. 
The corresponding results for the total cross section are presented in Table [4] and for the 
transverse momentum and rapidity distributions in Fig. [9] As we can see, the total cross 
section and the rapidity distributions from KrkNLO agree very well with those of MCFM 
and MCONLO. The difference between the KrkNLO results and the pure fixed-order ones is 
below 1% and comes from the 0(a s 2 ) contamination due to the MC PDFs, c.f. Eq.(4.4). 
Also, the differences between the a s (q 2 ) and a s (M|) are at the per-mille level. 

On the other hand, there are remarkable differences in the px,z distributions, especially 
at the lower end. They compensate in the total cross section between the first and the 
following bins. However, the choice of the running or non-running a s is numerically 
important for the px distribution at small and moderate values. The perfect agreement 
of MCONLO and MCFM at higher px is of course enforced by construction. This is not 
the case in KrkNLO, where some part of higher-order effects, beyond NLO, is included 
(which is also the case in POWHEG). We see, however, that both MCONLO and KrkNLO 
coincide below ~ 80 GeV in the case with a s (q 2 ). 

In the results of Fig. [9] the parton shower was artificially stopped after the first emis¬ 
sion, in order to make the meaningful comparisons with MCFM and limit the effects of the 
recoil due to subsequent emissions. In Fig.[l0]we lift this limitation and the parton shower 
goes to the very end, as in the normal operational mode of PSMC. As we see, switching 
on to the full PS influences considerably the low px part of the spectrum, in spite of a 
negligible effect on the rapidity distribution (and hence, on the total cross section). 

Again, KrkNLO with the running a s agrees very well with MCONLO at low and mod¬ 
erate px-, which results from the domination of the parton shower contribution in that 
region. However, the two KrkNLO results stay above MCONLO at higher px, which we 
again attribute to the admixture of some NNLO terms in the former. The difference 
between KrkNLO a s (q 2 ) and KrkNLO cr s (M|) at high px comes purely from the running 
of the coupling. 


5.2 All channels 

In the second round of the numerical tests, all the initial-state parton combinations con¬ 
tributing at NLO, that is the qq and qg channels are included. The total cross sections 
from KrkNLO, MCFM and MCONLO are compared in Table [5} As we can see, the differ¬ 
ence between MCFM/MCONLO and KrkNLO, which comes from the partial inclusion of 
the higher-order effects in the latter, is about twice as big as in the qq channel only. As 


discussed in Section [4721 this is related to the large gluon luminosity which leads to sizable 
differences between the MC and MS PDFs, and part of that difference shows up in the 
total cross section. The results in Table [5] are generally consistent with the differences 
staying below 5%. 


In Fig. 11 we show the distributions of the transverse momentum and rapidity of the 
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8 TeV : qq channel (full parton shower) 


8 TeV: qq channel (full parton shower) 




Pt,z [GeV] y z 


Figure 10: Comparisons of the transverse momentum and rapidity distributions from MCONLO 
and two versions of KrkNLO for the qq channel only as in Fig. [9] but the parton shower 
backward evolution runs to the end as normally. 


8 TeV: qq and qg channels (1st emission only) 




-3-2-10123 


Pt,z [GeV] 


Vz 


Figure 11: Comparisons of the transverse momentum and rapidity distributions from MCFM, 
MCQNLO and two versions of KrkNLO, now for both the qq and qg channels, with a s set 
as in the previous figures and the parton shower backward evolution stopped after the first 
emission as in Fig. [9] 
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[pb] 

MCFM 

1086.5 ±0.1 

MCONLO 

1086.5 ±0.1 

POWHEG 

1084.2 ±0.6 

KrkNLO a s (q 2 ) 

1045.4 ±0.1 

KrkNLO a s (M 2 z ) 

1039.0 ±0.1 


Table 5: Values of the total cross section with statistical errors for the Drell-Yan process, 
both channels, from the KrkNLO method compared to the fixed-order result of MCFM and 
the results of POWHEG and MC@NLO. 


Z-boson with the parton shower stopped after the first emission. The lower panel contains 
the ratio with respect to MCFM. Similarly to the case of the pure qq channel, we see the 
expected agreement between MCONLO and MCFM at large pr,z- The KrkNLO results 
with the fixed a s are not far from the MCFM ones in that region as well. On the other 
hand, all the PS-matched results differ from those of MCFM at smaller px, because of the 
lack of the Sudakov resummation in MCFM, see below for more discussion on this issue. 
Each curve behaves slightly differently but the variations are moderate. The rapidity 
distributions shown in Fig. 11 (right) are close to each other for all calculations. The 
~ 5% difference between KrkNLO and MCFM is the same as for the total cross sections 
in Table [5j In Fig. [Tl] we have also included the scale-uncertainty bands, which were 
produced using MCFM by means of changing (independently) both the factorization and 
the renormalization scales by the customary factors of 2 and 1/2. The differences between 
MCFM and KrkNLO are enveloped by these bands. 

Fig. 12 shows the similar distributions but, this time, the parton shower is allowed for 
an arbitrary number of emissions. 

There, we also show the results from POWHEG, as implemented in Herwig++. Even 
though the evolution variable is different (kx rather than q of Eq. (2.7) in all the other 
PS-matched results), POWHEG agrees quite well with the KrkNLO a s (q 2 ) version. This 
shows that the choice of the evolution variable in PS between kx and q is not so important 
for numerical results. The distributions of the transverse momentum and rapidity look 
quite similar to the previous results for the pure qq channel in Fig. 10 The transverse 


momentum distribution from KrkNLO agrees with MCONLO at high p TZ (up to a constant 
factor in the case with the fixed a s ), while larger differences are seen in the low px region. 

The differences of the order of 10-20% seen at low px between MCFM, MCQNLO, 
POWHEG and KrkNLO in Figs. 10, 11 and 12 reflect mainly the parton-shower feature of 
the soft-gluon resummation in the case of the matched results and the lack of it in the 
case of the fixed-order NLO calculations. More precisely, parton shower (without intrinsic 
kx ) has a gap below the minimum px ending the backward evolution and a non-physical 
spike ~ 3(px)- Normalization of the spike is governed by unitarization enforced by the 
shower. Kinematics, histogramming, intrinsic kx, etc. are smearing the shape of this non- 
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Figure 12: Comparisons of the transverse momentum distribution from MCQNLO, POWHEG 
and two versions of KrkNLO for both the qq and qg channels as in Fig. [TT] but the parton 
shower backward evolution runs to the end as normally. 


physical structure. Different implementations of the argument of a s strongly influence 
this structure near px = 0 and may show up in the ratios of the px distribution as large 
effects there. In particular, our results with cc s (M§) at low px are very different from 
the results with a s (q 2 ). The latter choice better represents resummation of subclasses of 
higher-order QCD contributions, relevant at low px, and it is therefore commonly adopted 
in PS MCs. In Fig. 12 another curve is added for MCQNLO with cc s (M§), demonstrating 
that the same kind of spike and dip near px = 0 is present in the corresponding ratios, 
even for the MCQNLO alone. Hence we conclude that the numerical results from KrkNLO 
near px = 0 are well understood and look as expected. 

Even better agreement is found in the comparison of results from KrkNLO a s (q 2 ) and 
POWHEG. Here, the two px,z distributions almost coincide in the whole range shown in 
Fig. 12 (left bottom panel). This is related to a similar, multiplicative, way of applying 
the NLO correction to the parton shower in KrkNLO and POWHEG. In both cases, this 
corresponds to the virtual part of the NLO correction being applied over the entire px,z 
range. This is different in MCONLO, where the virtual correction is spread by the shower 
only up to the upper limit of shower’s evolution variable (starting scale for BEV shower). 
That upper limit is typically set at the scale of a vector boson mass. Hence, in MCQNLO, 
the px,z distribution in the region px,z > mz is corrected only with the real part of 
NLO, and therefore it recovers exactly the fixed order NLO result in this region. The 
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Figure 13: Comparison of the rapidity distribution from MCFM, MCQNLO and two versions 
of KrkNLO for qq and qg channels as in Fig. [Tl] but parton shower backward evolution runs 
to the end as normally. 


multiplicative nature of the KrkNLO and POWHEG methods leads to mixed real-virtual 
corrections of the order a s 2 at high pr,z- These corrections are however part of NNLO, 
hence they go beyond the accuracy of NLO+PS result. 

The rapidity distributions shown in Fig. 13 agree very well in the central region be¬ 
tween all the PS-matched results, up to the normalization differences, the same as in 
Table [5] Differences at forward and backward rapidities come from increasing differences 
between PDFs in MC and MS schemes as x —» 1, c.f. Fig. [t} 

Finally, in Fig. 14 (left) we compare the KrkNLO results for the Z-boson transverse 
momentum distribution, from both channels and with full PS, with the corresponding 
result obtained from the fixed-order NNLO calculations implemented in the DYNNLO 
The NNLO distributions were obtained with the NNLO MSTW2008 PDFs 
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program 

and using pp = Pr = rriz, varied by the customary factors of 1/2 and 2 to estimate the 
uncertainty from neglected higher orders. Similar comparisons between the NNLO results 
and those of POWHEG and MC@NLO are shown in Fig. 14 (right). All the results are 
divided by the NLO distribution from MCFM. 

The calculations within the KrkNLO method were performed with two choices of the 
argument of the strong coupling a s : q 2 and M§, as was done earlier in the paper. These 
results are shown in Fig. 14 (left) as a solid blue line and dashed green line, respectively. 
In addition, in solid orange, we show the result of the KrkNLO matching in the case where 
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Pt.z [GeV] Pt.z [GeV] 


Figure 14: The Z-boson transverse-momentum distributions from KrkNLO compared with 


the fixed-order NNLO result from the DYNNLO program [431 (left). Similar comparisons for 
POWHEG and MCQNLO are also shown (right). All distributions are divided by the NLO 
results from MCFM. 


the argument of a s is set to the minimum of q 2 and M§. This is a physically motivated 
choice as at high p TZ our matched result should be driven by the fixed-order component, 
whose natural scale is M z . 


As we see in Fig. 14 (left), both the KrkNLO and the NNLO results show the same 


trends, quickly raising above the NLO result at low and moderate pr,z and then staying 
above it at high px,z (the curve for the case with a s (q 2 ) falls down due to the fact that 
the strong coupling is running even at high Pt,z)- The fact that our method gives the 
result that is higher than the NLO one at high px,z is a consequence of the mixed real- 
virtual 0(a s 2 ) terms, which constitute part of the NNLO correction and arise because of 
the multiplicative nature of the KrkNLO approach. Especially the a s (mm(q 2 ,m^)) choice 
does a good job, reproducing very closely the full NNLO distribution. Therefore this 
argument of the strong coupling is likely to be adopted in the future phenomenological 
studies. 


In Fig. 14 (right) we show similar comparisons with NNLO for MCONLO and POWHEG. 


The behavior at low px,z is close to that from KrkNLO. At high pr,z, however, MCONLO 
and POWHEG converge by construction to the NLO results, departing from the NNLO 
predictions. 


6 Summary and outlook 

We have discussed the KrkNLO method of matching the LO parton shower with the fixed- 
order NLO QCD corrections. The method is based on two elements: the change of the 
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factorization scheme from MS to the MC scheme, and upgrading the hardest emission 
to the full NLO accuracy by reweighting with a simple, positive weight. Details of the 
method and, in particular, demonstration of its NLO accuracy have been elaborated on 
in Section [3] 

The change of the factorization scheme allows one to eliminate troublesome ^-dependent 
terms from the coefficient function and, effectively, it amounts to creating the MC PDFs. 
In Section |4[ we have discussed how such PDFs can be obtained and how they differ from 
the standard MS parton distributions. There, we have also validated the MC factorization 
scheme by studying the Drcll-Yan process at the fixed-order NLO level and showing that 
the MS and MC scheme results are identical up to the order 0(a s ). 

We have implemented the KrkNLO method on top of the Catani-Seymour type of the 
parton shower in the Sherpa event generator for the case of the Z/ 7 *-boson production 
process (hence, the initial-state parton shower). In Section [5j we have presented the 
comparisons of the NLO-PS matched results obtained with our technique with the fixed- 
order NLO results from MCFM and with other matched results, namely those of MCONLO 
and POWHEG. In particular, we have demonstrated that the KrkNLO results recover the 
fixed-order NLO predictions (up to sub-percent differences for the qq channel only and 
~ 5% for all channels, coming from the beyond-NLO terms). 

As for the comparisons of KrkNLO with MCQNLO and POWHEG at the level of dif¬ 
ferential distributions, all three methods turn out to give essentially identical results for 
the yz spectrum. The px,z distributions look somewhat different for each method and 
the exact features depend on the initial channels and the recoil schemes. In general, the 
KrkNLO method provides similar predictions to the other two well-established approaches. 
In particular, the results with both channels and the full parton shower stay very close to 
those from the POWHEG method implemented in Herwig++. Residual differences come 
from spurious 0(a s 2 ) terms, which are different in each of the three matching methods. 

Applying the KrkNLO method to the Higgs production process is now under devel¬ 
opment and will be reported in the next paper. Applying the same KrkNLO method to 
more processes is quite straightforward, in the sense that the modification of the PDFs 
(to the MC scheme) and subtracting the hard process with the MC scheme counter-term 
(instead of the MS counter-term) can be done in the usual way for any process, also 
with more colored particles in the final state. The interesting question is rather whether 
this method eliminates all the non-physical ~ 5(pt) singularities, which are incompatible 
with the physical phase space of any PS MC? Following the study of Ref. |l5j, we are 
confident that for color-singlet heavy object production in the s and t channels this is 
true, although a general formal proof would be welcome. For more colored particles in 
the final state, this is still an open question. However, it is quite likely that the KrkNLO 
method reduces significantly the number of such pathological terms, and therefore will 
facilitate matching parton shower with any NLO and/or NNLO corrected hard processes. 
This is why, in our opinion, it is definitely worth to pursue this new development path in 
matching techniques. 

The current methods of the type NNLO+PS 


16-18*42 represent a clear progress in 


matching of the fixed-order NNLO QCD calculations with PS MCs, but still suffer from 
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various limitations; for instance, they are limited to a certain class of observables only. 
We are definitely thinking about extending KrkNLO to NNLO+NLOPS, in which NLOPS 
is a parton shower MC implementing the NLO evolution kernels in the fully exclusive 
form, thus providing the full set of the soft-colliuear counter-terms for the hard process. 
Ref. 15 reviews several feasibility studies, which show that constructing such a NLOPS 


is, in principle, plausible. In our opinion, any simplifications of the NLO+PS matching, 
as in the KrkNLO method, will be instrumental and very useful towards more ambitious 
fully exclusive NNLO+NLOPS projects. 
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APPENDIX 


A First emission in backward evolution 


In case of two LO parton showers in the BEV algorithm the cross section with exactly 
zero gluons in both showers includes two non-emission ^-functions 


(T LO - 

°o — 


dx F dx B dn ^(s,e) e- A *c(* F \°rf)-A&c(*B\ D^ c (s\s,x f ) D^ c (s\s,x b ), 


(A.l) 

where s = sx F x B . Using the identity of Eq. (3.30) twice, the above transforms immedi¬ 
ately into 




LO 


da 


= / dx F dx B dO — (s,6) e 5mc( * |s ’ 9s) SMc{S{s,qs) D^ c (s\q 2 , x F ), D^ c (s\q 2 , x B ), (A.2) 


dfl 


which coincides exactly with the n — 0 result in the forward evolution, starting from 
q 2 = q 2 B = q 2 . NB. In the case n = 0, we identify x F = x F and x B = x B . 

In the following we shall simplify the notation omitting the s argument in -Dmc(<? 2 > x) = 
DMc(s\q 2 ,x) and Smc((? 2 > Qi) = *S , Mc(s|q , |, q 2 ), wherever it is unambiguous. 

The distributions from the BEV algorithm with at least one gluon (n F +n B = 1, 2, 3...) 
reads as follows: 



where z\ integration limits are imposed by s(l — Z \) 2 / Z\ > q\ in the K kernels. 

The “unitarity” sum rule a^+a^ = a LO is of course automatic in the BEV algorithm. 
Nevertheless, proving it provides an interesting cross-check. The differentiation over the 
lower boundary in the q 2 integral defining the Afunction leads to: 


d 


d In q 2 


D -A M c(x F \s,q 2 )-A MC (x B \s,q 2 ) _ 


dz 


= / ^[K MC (x F \z,q 2 ) + K MC (x B \ z ,q 2 )] e -^McU F b, 9 2 )-^Mc(i B l^ 2 ). 


(A.4) 
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After combining two integrals in Eq. (A.3) we obtain 


In s 


= J dx F dx B di1 J dln(q , 1 ) 
i n?f 

da 


d 


d In q\ 


g l s >9l) l s '9l) 


X dQ^^ 3 ' 9 ^ ° MC ( S,£f ) ^(s-^) 


= / dx F dx B dQ 


X _ e - 4 Mc(Al s * 2 )"^Mc( i fll s ’ , ! s 2 ) 


(A.5) 


from which the sum rule <jg ° + a^+ = a LO results immediately. 

Before we transform Eq. (A.3) into the FEV picture, let us first expand K-kernel^J 


a\2 — / dx F dx B dfl 


*< ifjf 

■ „2 F T 

9s ■*? 


S _v 2 1 


£7 ( q 2 “F 

2 c -^mc(^fI s .9i f ) c -^mo( £ bI^ 9? f ) MC V 1f ’ Z1 

S (l- Zl )7zi>g2 e Fe 7VF , 2 A 

w mcWi F )A 




7 gf fl d *1 

9? *8 


da 


x Dm C (s,x f ) A>mc(s,£ b ) —(s£ f £ b ,0). 

Now using again the identities 

^mc(^fI s, 9i f ) _ r:( 5 »9?j.) ^MC (?lf ? *^f) — ^mc(*b I s >9i f ) __ r .~£Mn( s i9T F ) ^-b) 

^Mc(hA) 

we arrive at the following FEV representation 

da 
df2 


(A.6) 


^Mc( s ) %f) 


a±° = / dXpdx B df2 ——(sx F x B ,§) 


s 1 

„2 


^ fdZi 

* ' ' „2 I Vs(l-2i) 2 /2l>9 2 


Tip J Zl 

■ „2 T 

9s X F 


Td^ c {<,1,x b ) 

F Zl 


1 


dq\ B [ dz l 


{ Ql B J z i 

<& x b 


— ds(l-z 1 y/z 1 >g 2 1 


T(A)e s mc(^ b ) sg c ^ D ^ c { q \ x F ) D ^ c ( q \J^-) 

1 

(A.7) 


We denote T(z) = 


TV 1 — 2 
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Finally, the longitudinal integrations are streamlined with the help of the substitutions 
to the x -variables before the emission, x F /z\ —> x F , X\ B = x B or x B fz\ —> x B , X\ F = x B1 
respectively in each shower, and the convolution structure is made manifestly symmetric: 


— I dx I dx F dx B df2 


x 


S 1 

f dq\ f dz l 




^s(i- 2 i) 2 / 2 i>gj 'P(^i) e MC ^ ’ 91f ' ) mc ^ ,qiF ' > D^ /[C (q lF , x F )D^ /[C (q lF , x B ) 

Z\ F 


"F 

1 


dqf f dzi 




3>(zi) e 


S Mc( S ’l1 B ) , S'mc( S ’ i J1h) F) f (n 2 ry ^ DB (rd 


1b -^Mc(? 1 r’ X f)Dmc(, C Li b i X b) 


X x ~ Z l‘ X F X B) 


(A.8) 


The above formula is for the FEV algorithm but for one gluon next to the hard process and 
any number of trailing gluons down to q 2 s . The same formula in Eq. (3.33) was obtained 
by means of factorizing such a gluon from the fully exclusive FEV formula of Eq. (3.24). 


B Exclusive NLO corrections in MS scheme 


Following Ref. |8], the bare NLO real + virtual unintegrated DY cross section for the 
incoming qq with the effective mass of ybq, in d — 4 + 2e dimensions, is quite simple: 


o%£{s i,e)= I dQ J dadfl {S(a)<S(/3) W 0 ( Sl ) + W 2 ( Sl , a, fd) 

a+/5< 1 

W 0 (si) = (l + V(si,e))-^y(si,6 l ), 

W 2 (si, a, f3) = H 2 (si,a,/3,e) A) + H 2 (s 1 ,/3,a,e ) j^{ zs i A) 


V(si,e) = 


2Cfol s (Att) e r{l + e) ( S\\ e { 1 3 1 7 t 2 \ 


i7 2 (si,a,/3,e) = 


7T r(l + 2e) Vh 2 J l 2e 2 + 4e 2 12 J ’ 

2Cfoi s (47r) _e f sia(3 V (1 — a) 2 + ea 2 


it r( 1 + e) \ /i 2 / 


2 a/3 


(B.l) 


where z = l — a — /3 and notation of Section 3.1 is used. The integration over the azimuthal 
angle is already done in order to improve readability. The integration decay lepton angles 
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and a and /3 keeping z = 1 — a — (3 fixed, leads to a divergent (bare) result 


cr 


Sb( s i» c )= / dza 0 (siz)[5(l-z)+p B (z,s l ,e)\, 


'0 


2 a s (4tt 

Pb{z,s i,e) = - 

7T 




1 + £ 2 


-2J 5(1- 
+ (1 + * 2 ) 


(B.2) 


2 , ln(l - *) 


,2(1 — z) J + ' - ' " 7 V 1 
The MS recipe tells us to subtract from the bare radiator function p B its pole part 


(e) = PP p B (z, Sl ,e) = —C F - 


2a., 


1 (47t)~ 


MS 


7r 


e T(1 + e) \2(1 — z) 


1 + z 2 


(B.3) 


After including PDFs and integrating over Si which leads to the replacement Si —> sx F x B , 
we obtain the finite hadron-level cross section 


a 


NLO 


(s) = j dx F dx B dz |h(l — z) 

+ [(1 - PP )p B (z,sx F x B ,e 


(B.4) 


szx F x B )D^ b (p~, x F )Df b (p^ x B ). 


The contribution from the qq channel to the inclusiv e coef ficient function of the Drcll-Yan 

The above formula can be used 


process is 
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3.11 


= (l-PP )p B (z, of Eq. 

to define The PDFs and dig them out from experimental DY data. Alternatively, it can 
be used to predict the DY cross section taking the PDFs obtained from other processes, 
e.g. deep inelastic lepton-proton scattering (D1S). 

The above formulation is suitable for the total cross section. In the case of any NLO- 


class exclusive observable defined by the function J = Jnlo{x f , x b , z, kf T ) of Eq. (3.35), 


the above formulas are reorganized using the technique of the Catani-Seymour soft- 
collinear counter-terms (SCC). In the present case SCC coincides (for e = 0) with the 
single gluon LO distribution of the parton shower MC: 

2 C F a s (47r) _e f Sia/3 V1 + z 2 + e(l — z) 2 da 0 


W 0t ( Sl ,a,P,e) = — WT7) [-r) 


it l (i + e) \ p* j 2a/3 

where z = 1 — a — (3. The above SCC is subtracted and added: 


dn (zs " e) ’ (R5) 


cr 


NLO 


(s)[j] = / dx F dx B Df s (p 2 ,x F )Df s (p 2 ,x B ) 


x (l-PP) 


J(x F ,x B , 1,0) 


dft W 0 (si) + J dad/3 J dft W c t(si, a, /3, e) 

Oi-\-(3< 1 

dfl Jdad/3 [W 2 (si, a, (3)J(x F ,x B ,z, k\ T ) - W rf (si, a, f$) J(x F , x B , 1, 0)] £=0 }>, 

a+/3<l 


(B.6) 


33 The assignment p 2 = zs\ induces the presence of the term ~ ln(z)(l + z 2 )/{ 1 — z) in C%l s . 
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where Si = sx F x B . In the second real emission integral W ct eliminates completely the 
soft/collinear singularities, so it can be evaluated in the e —» 0 limit. The key point 
here is that W ct is multiplied by Jlo = J(x F , x B , 1,0), hence the partial integration 
J dad/35(z — a — (3) can be performed for the fixed z in order to isolate and subtract the 
1/e pole: 


cr NLO (s)[J] = J dx F dx B dft J dz D^ is (/j 2 ,x f )D^ is (/i 2 ,x b ) 

x {( 1 + A vs ) <5(f “ z ) +2Z q (z)} C ^{zs 1 ,6)J(x F ,x B ,l,0) 

+ J dx F dx B dfI D^ s (fi 2 , x f)D^ s (i^i 2 , x b) J dad/3 AW rea i[J], 

ct+/3<l 


with 


2 Z g (z) = 


2 Cpa s 


7r 


1 + z 2 


In 


(l-,) s 


+ 


1 + z 2 


In — + 


1-z 


2(1 — z) z 2(1 -z) /i 2 ' 2 

AW real [J] = [W 2 (si, a, (3) J(x p , x g , k\ T ) - W ct (si, a, (3) J(x F , x B , 1,0)] £=0 , 


(B.8) 


where A ( y S = 2C ^ as (^- — |) is that of Eq. (3.10), U q (z)\s=^ = AC 2q , see Eq. (3.13), and 
s = szx F x B as usual. 

In the MC scheme the entire ~ 5(k'l T )S(z) part gets eliminated (modulo 0(a 2 ) terms) 
thanks to the assignment s = /i 2 and redefinition of the PDFs 


D“ ? c (/i ,:r) = / dzdx'5(x - zx')[5( 1 - z) + S q (z )\ s=fl 2 (B.9) 


see Eq. (3.17). The only freedorcpjis in the coefficient of 5(1 — z). However, the convention 
f dzE(z) = 0 adopted in Ref. |19| re moves this freedom and fixes the content of Ays- 
Inserting Eq. (B.9) into Eq. (B.7), we obtain the formula for the NLO cross section in 
the MC scheme 


a 


NLO 


(s)[J] = / dx F dx B di1 j dz D^ ic (s,xf)D^ [C (s,xb) 


} mc, 


X { 1 + A Vs)j^( S ^ 9 ) J ( X f, x bA,0) 


(B.10) 


dx F dx B dQ I dad/3 AW rea i[J] D^ ic (s, xf)D^ ic (s,xb), 


'.MC/ 


q+/5< 1 


where we have replaced the MS PDFs by the MC PDFs also in the real part. The cross 
sections in Eqs. (B.10) and (B.7) are equal up to 0(a 2 ). 


34 The z + 0 part of S is unique it has to be eliminate completely from the NLO distribution, because 
the singular non-positive ~ 5(k\ T )S term cannot be included in the multiplicative MC weight. 
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Finally, let us remark, following Ref. 


19 


that formally the transition from the MS 


to MC scheme can be done by means of subtracting from the bare NLO distribution the 
term 


ditto = ®f) + z,(z)\t=,*, 


_ pi 1 ] 


(B-ll) 


instead of the pure pole /^(e). It is appealing, that /^(e) concides with the integrated 
single-gluon emission distribution of PS MC extrapolated to d = 4 + 2e dimensions. This 
justifies the name of the MC factorization scheme. 
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